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Introduction 


Despite  advances  in  the  technologies  for  the  diagnosis  and  treatment  of  breast  cancer,  strategies 
for  the  prevention  of  this  disease  are  poorly  defined.  However,  one  fundamental  concept  arising 
from  recent  studies  of  mammary  gland  biology  is  the  understanding  that  the  processes  involved 
in  normal  mammary  development  and  carcinogenesis  are  intrinsically  related  (1).  One  of  the 
most  frequently  cited  examples  of  this  principle  is  the  protective  effect  of  an  early  full-term 
pregnancy  against  breast  cancer.  There  is  strong  epidemiological  evidence  that  women  who 
experience  a  full-term  pregnancy  early  in  their  lives  have  a  significantly  reduced  risk  of 
developing  breast  cancer  (2,  3).  This  is  substantiated  by  studies  in  a  rat  model  in  which  prior 
treatment  with  estrogen  and  progesterone,  (to  simulate  the  effects  of  pregnancy)  confers 
resistance  to  tumorigenesis  in  animals  subsequently  challenged  with  chemical  carcinogen  (4-12). 
The  E+P-treated  rat  is  thus  a  well-established  model  to  study  the  effects  of  parity-induced 
protection  against  mammary  cancer.  In  spite  of  this  fact,  little  is  known  of  the  specific 
mechanisms  governing  pregnancy-specific  developmental  changes  in  the  mammary  gland.  In 
this  respect,  the  elucidation  of  novel  targets  for  estrogen  and  progesterone  action  in  the 
mammary  gland,  is  crucial  to  our  understanding  of  how  an  aspect  of  normal  development  might 
mediate  the  response  of  the  organ  to  future  proliferative  signals. 

To  address  this  question,  we  have  used  the  well-characterized  Wistar-Furth  rat  model 
(described  previously,  4),  in  conjunction  with  subtractive  suppressive  hybridization  (SSH) 
methodologies  (13)  to  identify  markers  that  are  persistently  up-regulated  in  response  to  prior 
exposure  to  E  and  P  in  the  mammary  gland.  As  a  result  of  these  initial  studies,  we  selected 
several  specific  markers  to  study  in  greater  detail  to  investigate  their  role  in  relation  to  parity- 
related  protection.  To  achieve  this  goal,  we  proposed  to  use  mRNA  expression  analysis  to  study 
the  temporal  and  spatial  expression  of  each  candidate  biomarker,  to  evaluate  their  role  as  a 
function  of  reproductive  history  in  relation  to  both  differentiation  and  parity -related  protection  of 
the  mammay  gland.  Furthermore,  since  the  role  of  many  of  these  markers  has  not  yet  been 
described  in  the  mammary  gland,  we  described  experiments  to  further  elucidate  their  function  in 
relation  to  parity-related  changes  within  the  gland.  In  part,  this  approach  will  involve  the  use  of 
novel  retroviral  transduction  methods  to  overexpress  the  gene  of  interest  in  a  parity-independent 
manner  in  the  reconstituted  mammary  glands  of  recipient  animals.  If  this  approach  is  succesful, 
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then  we  have  proposed  experiments  to  investigate  the  influence  of  overexpression  of  a  gene  of  interest 
in  relation  to  protection  against  tumorigenesis  following  carcinogen  challenge. 


Report  Body 

a)  Summary  of  research  objectives 

The  specific  objectives  of  this  project,  described  in  the  original  proposal  are  listed  below: 

Objective  1:  To  characterize  the  role  of  candidate  biomarkers  and  to  investigate  their  expression 
in  the  rat  mammary  gland: 

a)  To  determine  their  temporal  and  spatial  pattern  of  expression  in  the  rat  mammary  gland. 

b)  To  determine  whether  these  candidate  genes  are  markers  of  differentiation,  per  se  or  of 
protection. 

Objective  2:  To  investigate  the  function  of  candidate  genes  in  vivo: 

a)  To  investigate  the  function  of  GB7  and  to  determine  its  role  during  normal  cellular  processes 

b)  To  investigate  the  function  of  GB7,  RbAp  and  Napl  in  reconstituted  mammary  glands 

c)  To  determine  the  role  of  systemic  factors  on  parity-induced  changes  in  the  mammary  gland. 


Many  of  the  experiments  described  in  objective  1  of  the  proposal  have  been  completed,  and  the 
results  published  in  a  peer-reviewed  article  in  November,  2001  (14).  In  situ  hybridization 
analysis  was  used  to  determine  the  spatial  expression  of  G.B7  and  RbAp46.  The  same  approach 
was  also  used  to  determine  the  effect  of  perphenazine  treatment  on  the  expression  of  these  two 
markers.  A  quantitative  reverse-transcriptase  PCR  method  has  been  developed  to  examine  the 
temporal  expression  of  these  markers  in  the  rat  mammary  gland. 

The  experiments  described  in  Objective  2  are  in  progress.  However  as  you  may  be  aware,  the 
Tropical  Storm  that  struck  Houston  in  June  2001  caused  significant  impairment  to  researchers  at 
Baylor.  This  included  not  only  closure  of  the  research  facilities  in  the  immediate  aftermath  of 
the  disaster,  but  also  an  extended  period  of  downtime  due  to  the  losses  incurred  as  a  result  of  the 
disaster.  It  is  estimated  that  members  of  our  laboratory  lost  two  months  of  research  time  due  to 
the  constraints  imposed  due  to  lack  of  access  to  the  research  facility,  cataloguing,  ordering  and 
attempting  to  salvage  damaged  research  materials,  and  waiting  for  the  replacement  materials  to 
arrive  before  we  could  recommence  work. 
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b)  Results 


Objective  1 

a)  To  determine  their  temporal  and  spatial  pattern  of  expression  in  the  rat  mammary  gland. 

b)  To  determine  whether  these  candidate  genes  are  markers  of  differentiation,  per  se  or  of 
protection. 

In  the  original  statement  of  work,  we  proposed  to  use  Northern  analysis  and  in  situ  hybridization 
to  investigate  the  expression  of  three  markers  (G.B7,  RbAp46  and  Napl)  at  different  stages  of 
development  in  the  rat  mammary  gland.  In  addition,  we  proposed  to  use  perphenazine  treatment 
as  an  experimental  paradigm  to  distinguish  the  molecular  processes  involved  in  protection  of  the 
gland  from  those  that  are  inherent  to  pregnancy-related  differentiation  of  the  gland.  In  the 
previous  Annual  Report,  we  described  the  use  of  i n  situ  hybridization  analysis  to  determine  the 
spatial  expression  of  G.B7  and  RbAp46.  We  also  described  experiments  to  recapitulate  the  effect 
of  perphenazine  treatment  in  Wistar-Furth  rats  and  the  use  of  in  situ  hybridization  analysis  to 
determine  the  effect  of  perphenazine  treatment  on  the  expression  of  G.B7  and  RbAp46.  The 
results  of  this  anlysis  was  also  reported  in  a  recent  manuscript  (14).  We  also  described  the 
difficulties  of  detecting  low  abundance  mRNAs  such  as  G.B7  and  RbAp46,  the  necessity  of 
using  PolyA+  RNA  for  Northern  analysis  and  technical  difficulty  of  isolating  sufficient  amounts 
of  poly  A+  RNA  from  mammary  glands  at  different  stages  of  development  to  satisfy  this 
requirement. 

To  overcome  these  difficulties,  we  proposed  two  alternative  approaches  to  enable  us  to 
quantitate  the  expression  of  these  mRNAs  without  the  necessity  of  isolating  poly  A+  RNA  (as 
described  in  the  previous  progress  report).  In  the  first  case,  we  proposed  to  apply  the  use  of 
cDNA  microarray  technology  to  quantitate  the  expression  each  marker  at  different  stages  of 
development.  We  begun  experiments  to  investigate  the  effectiveness  of  this  approach  by 
constructing  a  pilot  array,  containing  96  non-redundant  clones  from  the  E+P  subtracted  SSH 
library  (including  G.B7,  Napl  and  RbAp46),  in  conjunction  with  the  Microarray  Facility  at 
Baylor  College  of  Medicine.  The  cDNA  inserts  of  these  clones  were  amplified  using  the  nested 
primers  in  the  SSH  library  backbone,  purified  by  gel  filtration  chromatography  and  arrayed  using 
silanized  glass  technology  by  means  of  the  384  well  format  array  printer  in  the  Microarray  Core 
Facility.  To  investigate  the  effectiveness  of  this  approach,  we  screened  the  pilot  array  with 
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fluorescently-labeled  first  strand  cDNA  probes  prepared  by  reverse  transcription  of  total  RNA 
isolated  from  mammary  glands  of  E+P  and  age-matched  virgin  (AMV)  rats.  We  conducted  these 
experiments  using  a  variety  of  hybridization  conditions  to  determine  the  optimal  conditions  for 
sensitivity  of  the  detection  method.  Unfortunately,  the  results  of  these  preliminary  experiments 
showed  that  less  that  30  %  of  the  cDNA  targets  (including  those  of  G.B7,  Napl  and  RbAp46) 
could  be  detected  by  this  method,  even  though  the  library  from  which  they  were  isolated  was 
constructed  from  the  same  E+P  RNA  pool  used  for  the  probe  synthesis.  We  believe  that  this  is 
once  again  due  to  the  inherent  difficulties  of  detecting  low  abundance  RNAs  from  a  relatively 
quiescent  tissue  such  as  the  involuted  mammary  gland.  Accordingly,  we  began  to  investigate 
alternative  methods  of  probe  preparation  using  indirect  labeling  with  Cy3/Cy5  or  cDNA 
amplification  to  increase  the  sensitivity  of  detection.  However,  in  the  last  year  the  Microarray 
facility  at  Baylor  ceased  its  array  screening  services  to  reevaluate  the  method  of  array 
production. 

As  an  alternative  approach,  we  decided  to  apply  the  technique  of  real-time  PCR  analysis  to 
investigate  the  expression  of  these  markers  at  different  stages  of  development,  and  as  a 
consequence  of  protection  or  differentiation  of  the  gland.  This  approach  has  the  advantage  of 
being  extremely  rapid  and  sensitive  to  small  changes  in  the  concentration  of  a  reverse- 
transcriptase  amplified  cDNA  product.  The  application  of  this  technique  was  previously  limited 
by  lack  of  access  to  the  facilities  for  performing  this  analysis  and  the  cost  of  purchasing 
commercial  reagents  such  as  the  double-labeled  fluorescent  probes  and  proprietary  amplification 
kits  required  for  these  experiments.  However,  recently  our  department  has  acquired  the  facilities 
for  performing  quantitative  real-time  RT-PCR.  Furthermore,  a  number  of  recent  publications 
have  described  the  use  of  SYBR  Green  (as  opposed  to  ‘Taqman’  or  FRET-based)  detection 
methods,  combined  with  self-assembled  PCR  reagents  (15, 16).  Thus  we  decided  to  adapt  this 
approach  to  study  the  expression  of  candidate  markers  in  the  rat  mammary  gland.  So  far  we 
have  applied  this  analysis  to  the  RbAp46,  and  another  potentially  interesting  marker,  Stral3  (17). 
Expression  of  Cytokeratin-8  was  used  as  a  normalizing  control,  to  correct  for  the  epithelial 
content  of  mammary  glands  at  different  stages  of  development.  Using  the  human  genomic 
database  sequences  (http://www.ncbi.nlm.nih.gov/genome)  and  a  search  algorithm  to  predict 
exons  by  aligning  genomic  and  cDNA  sequences  (18),  we  predicted  the  exon-intron  structure  of 
the  homologous  rat  genes,  and  then  designed  mRNA-specific  primer  pairs.  These  primers  either 
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spanned  an  exon-intron  boundary,  or  if  this  was  not  possible  then  the  two  primers  were  encoded 
by  separate  exons  of  the  gene  of  interest.  In  the  latter  case,  amplification  times  were  shortened 
to  prevent  the  amplification  of  a  long  genomic  fragment.  We  designed  the  primers  to  produce 
short  amplicons  (150-200  bp)  to  allow  rapid  cycling  conditions,  which  should  also  prevent  the 
amplification  of  non-specific  products.  Primers  were  then  rigorously  tested  by  BLAST  analysis 
to  confirm  the  specificity  of  the  primer  sequence  and  in  a  standard  RT-PCR  reaction  (using  no- 
RT  and  genomic  DNA  as  controls  and  the  same  amplification  conditions  as  the  quantitative 
assay)  to  ensure  the  absence  of  any  contaminating  genomic  products.  The  resulting  amplicons 
were  analyzed  by  gel  electrophoresis  on  2  %  agarose  gels  to  confirm  specificity  of  the  reaction 
and  the  absence  of  primer-dimers.  Synthetic  RNA  standards  were  also  incorporated  into  the 
assay  to  allow  absolute  quantitative  of  the  levels  of  the  gene  of  interest  (19).  Synthetic  standards 
were  generated  by  in  vitro  transcription  from  a  vector  possessing  a  T7  promoter  adjacent  to  the 
sense  strand  of  the  cDNA  target.  The  transcript  was  quantitated  by  incorporating  trace  amounts 
of  [a32P]  rUTP  into  the  transcription  reaction  and  the  transcript  analyzed  by  polyacrylamide  gel 
electrophoresis  to  ensure  that  full-length  transcripts  are  produced.  The  use  of  an  RNA  standard 
instead  of  a  DNA  standard  circumvents  any  problems  due  to  differences  in  amplification  as  a 
result  of  the  reverse-transcription  step. 

The  assay  was  then  conducted  as  a  two-step  process:  the  first  step  involving  reverse- 
transcription  of  RNA  in  a  standard  thermocycler,  and  in  the  second  step  the  reverse-transcription 
products  were  simultaneously  amplified  and  the  amplification  products  detected  using  an  ABI 
Prism  7700  Sequence  Detection  System.  Duplicate  samples  were  analyzed  in  same  assay  along 
with  serial  dilutions  of  reverse-transcribed  synthetic  standard  (as  described  above).  At  the  end  of 
the  reaction,  amplification  plots  were  inspected  to  make  sure  that  none  of  the  templates  had  been 
overloaded,  that  a  single  amplification  product  was  produced,  and  that  replicate  samples 
produced  overlapping  amplification  curves  (figure  1).  A  standard  curve  was  generated  by 
plotting  Cj  values  for  each  concentration  of  standard,  versus  the  number  of  target  molecules  in 
that  reaction.  The  linearity  of  the  response  was  examined  by  calculating  the  regression 
coefficient  for  the  plot  (figure  2).  As  stated  above,  we  have  used  this  approach  to  examine  the 
expression  of  RbAp46  and  Stral3  at  different  stages  of  development  and  as  a  consequence  of 
differentiation  or  protection  of  the  gland.  We  hope  to  extend  these  studies  to  G.B7,  using  the 
recently  assembled  draft  mouse  genomic  sequence  database 
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(http://www.ensembl.org/Mus_musculus)  to  predict  the  exon-intron  stucture  of  this  gene  in 
mouse  and  rat. 


Figure  1.  Amplification  plot  for  varying  concentrations  of  RbAp46  standard 


-2.854x  +  34.741 


R2=0.995 


In  the  final  part  of  this  objective,  we  proposed  to  continue  to  characterize  candidate 
markers  by  sequencing  and  Northern  blot  analysis.  In  addition  to  the  data  published  in  our 
manuscript  in  November,  we  have  also  sequenced  96  clones  from  a  reverse  subtracted  (E+P 
driver,  AMV  tester)  cDNA  library  which  was  constructed  and  screened  using  the  same  approach 
described  for  the  E+P  subtracted  library. 


Objective  2 

a)  To  investigate  the  function  of  GB7  and  to  determine  its  role  during  normal  cellular 

processes. 

In  the  first  part  of  this  objective,  we  proposed  experiments  to  determine  the  function  of 
GB7  and  to  predict  its  role  in  normal  cellular  processes.  As  started  in  the  original  proposal  and 
previous  progress  report,  GB7  is  a  novel  cDNA  whose  pattern  of  expression  suggests  it  is  an 
excellent  marker  for  parity-induced  protection.  In  spite  of  this  fact,  we  have  been  unable  assign 
a  function  to  this  marker.  In  the  previous  progress  report,  we  described  the  isolation  and 
sequencing  of  full-length  clones  for  this  mRNA.  We  then  conducted  extensive  database  analysis 
in  an  attempt  to  elucidate  the  function  of  G.B7;  however,  this  search  yielded  little  further 
information  other  than  the  identification  of  a  potential  human  homologue  for  this  gene,  which 
mapped  to  a  region  spanning  a  known  chromosomal  break-point  associated  with  a  number  of 
human  cancers  (20). 

Recently,  we  decided  to  repeat  this  analysis  with  the  hope  if  identifying  newly 
characterized  human  ESTs  with  similarity  to  G.B7.  The  full-length  G.B7  mRNA  sequence 
(Genbank  Accession  #  AY035343)  was  subjected  to  lower  stringency  alignment  computation 
using  the  analysis  algorithms  of  Pipmaker  and  Multipipmaker.  In  the  first  case,  the  full-length 
rat  mRNA  sequence  was  aligned  with  the  sequence  of  the  entire  human  BAC  pertaining  to  the 
region  of  homology.  In  the  second  analysis,  sequence  from  two  overlapping  rat  BACs  with 
homology  to  G.B7  were  added  to  the  alignment.  These  additional  sequences  were  added  to 
strengthen  possibility  of  producing  a  good  alignment,  since  the  first  alignment  required 
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homology  between  orthologous  mRNA  and  genomic  sequences.  This  analysis  was  performed  in 
collaboration  with  Dr.  Monique  Rijnkels,  a  senior  postdoctoral  fellow  in  Dr.  Rosen’s  laboratory. 
The  results  of  this  analysis  confirm  the  original  prediction  of  a  homologous  human  G.B7 
sequence,  located  on  human  chromosome  2.  Furthermore,  it  suggested  that  the  putative  human 
gene  is  encoded  for  the  most  part  by  a  single  large  exon  of  4503  bp  and  two  smaller  exons  of  147 
and  376  bp  respectively.  This  information  should  enable  to  predict  suitable  regions  to  design 
primers  for  real-time  PCR  analysis  (as  described  above).  Finally,  the  results  of  this  analysis 
showed  that  a)  the  G.B7  coding  region  overlaps  with  one  exon  of  the  ORF  of  a  hypothetical 
human  gene;  and  b)  that  it  is  close  to  a  predicted  imprinted  region  (22).  As  we  have  stated 
previously,  we  believe  that  this  gene  may  function  as  a  non-translated  regulatory  RNA  (14).  In 
the  last  few  years,  the  importance  of  non-coding  RNA  genes  has  grown  in  recognition,  with 
more  and  more  examples  of  this  class  of  molecules  being  discovered  (23).  Although  this  class  of 
molecules  may  be  of  diverse  function,  there  is  evidence  linking  a  number  of  non-coding  RNAs 
with  the  regulation  of  imprinting  of  certain  genes  or  loci.  The  fact  that  the  human  homologue  of 
GB7  overlaps  with  a  distinctly  separate  putative  OFR  suggests  that  G.B7  may  also  play  a  role  in 
tissue-  or  developmentally-specific  imprinting.  The  results  of  the  human  genomic  database 
analysis  may,  therefore,  enable  us  to  design  further  experiment  to  investigate  the  function  of 
G.B7. 

b)  To  investigate  the  function  of  GB7,  RbAp46  and  Napl  in  reconstituted  mammary  glands 

The  second  part  of  this  objective  proposed  the  use  of  reconstituted  mammary  gland 
experiments  to  elucidate  the  function  of  GB7,  RbAp  and  Napl  in  an  in  vivo  context.  To  achieve 
this  goal,  we  have  recently  developed  methods  for  growing  rat  mammary  epithelial  cells  in 
primary  culture  and  infecting  them  with  a  recombinant  retrovirus  expressing  a  gene  of  interest. 

In  these  preliminary  experiments,  we  have  tested  several  retroviral  constructs  and  one  lentiviral 
construct  expressing  a  p-galactosidase  reporter  gene  to  determine  the  efficiency  of  the 
transduction  with  different  vector  systems.  MEC  were  prepared  from  96  day  old  virgin  rats  and 
maintained  in  culture  as  previously  described  (24).  After  48  hours  of  culture,  primary  MEC 
were  infected  with  freshly  prepared  virus  stocks.  These  stocks  were  filtered  before  use  to 
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remove  residual  virus  particles  and  a  constant  volume  of  each  applied  to  48  h  primary  cultures. 
After  five  days  of  culture,  and  four  rounds  of  infection,  MEC  were  stained  with  the  chromogenic 
substrate  for,  X-gal,  and  counterstained  with  nuclear  fast  red  stain  to  estimate  infection 
efficiency.  Using  the  lentiviral  vector,  we  have  been  able  achieve  infection  of  five  to  ten  percent 
of  the  total  cell  population  —  measured  by  the  nuclear  localization  of  blue-staining  following 
incubation  of  the  cultured  MECs  with  X-gal  (figure  3).  This  is  five  to  ten  times  higher  than  the 
efficiency  of  infection  achieved  by  conventional  retroviral  strategies  employing,  as  an  example, 
the  Moloney  Murine  Leukemia  Virus  (MMLV)  derived  vector  pS2. 

In  subsequent  experiments  MECs  were  again  cultured  for  five  days,  and  during  that  time 
subjected  to  four  rounds  of  infection  with  the  same  (3-galactosidase-expressing  lentivirus  as 
described  above.  MECs  were  then  removed  from  the  culture  plates,  washed  several  times  in 
serum-free  media  and  injected  into  the  cleared  fat  pads  of  21  day  old  syngeneic  host  animals.  As 
a  control,  non-transfected  MECs  were  also  injected  into  the  contra  lateral  gland  of  two  of  the 
same  rats.  Eight  weeks  following  transplantation,  the  recipient  animals  were  sacrificed  and  the 
outgrowths  removed  and  stained  with  X-gal  to  detect  those  resulting  from  transduced  MECs.  In 
this  experiment,  two  out  six  outgrowths  from  the  transduced  MECs  resulted  in  blue-staining  out¬ 
growths,  indicating  that  we  can  produce  transduced  mammary  outgrowths  by  this  method. 
Future  experiments  will 
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Figure  3:  p-galactosidase  staining  of  primary  MEC  cultures  following  two  rounds  of  infection  with  recombinant 
lentivirus  vector  expressing  LacZ 
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Key  Research  Accomplishments 


•  Experiments  involving  perphenazine  treatment  have  been  reproduced  in  Wistar-Furth  rats. 
Level  of  morphological  development  examined  by  whole  mount  analysis,  in  comparison  to 
mammary  glands  for  E+P-treated,  pregnant  and  unstimulated  control  animals.  (Objective  1) 

•  In  situ  hybridization  analysis  used  to  characterize  the  expression  of  G.B7  and  RbAp46  in 
response  to  pregnancy,  hormonal  stimulation  and  perpehnazine  treatment.  (Objective  1). 

•  Screening  of  rat  E+P  SSH  library  complete.  Twenty-one  markers  from  this  library 
characterized  by  Northern  blot  analysis  to  verify  their  pattern  of  expression.  (Objective  1). 
Reverse  subtracted  library  constructed,  screened  and  96  clones  partially  sequenced. 

•  Manuscript  revised  and  accepted  for  publication  (July,  2001) 

•  A  quantitative  reverse-transcriptase  PCR  method  has  been  developed  to  examine  the 
expression  of  candidate  markers  in  the  rat  mammary  gland  (Objective  1). 

•  Isolation  and  sequencing  of  full-length  clones  for  G.B7  completed.  Sequence  submitted  to 
Genbank  (accession  #  AY035343)  (Objective  2). 

•  Identification  of  putative  human  homologue  for  G.B7  (Objective  2). 

•  Established  methods  for  culturing  mammary  epithelial  cells  (MEC)  from  rats  (Objective  2) 

•  Established  methods  for  transducing  rat  MECs  using  a  lentivirus  vector,  and  for  using  these 
transduced  MECs  to  produce  reconstituted  mammary  glands  expressing  a  p-galactosidase 
reporter  gene  (Objective  2). 


Reportable  Outcomes 


Manuscripts: 

Melanie  R.  Ginger,  Maria  F  Gonzalez-Rimbau,  Jason  P.  Gay  and  Jeffrey  M.  Rosen.  2001. 
Persistent  Changes  in  Gene  Expression  Induced  by  Estrogen  and  Progesterone  in  the  Rat 
Mammary  Gland.  Mol.  Endocrinol.,  Nov.  15  (1 1):  1993-2009  (Refer  to  Appendix  i) 

Medina  D,  Sivaraman  L,  Hilsenbeck  SG,  Conneely  O,  Ginger  M,  Rosen  J,  O’Malley  BW.  2000. 
Mechanisms  of  hormonal  prevention  of  breast  cancer.  Ann  N  Y  Acad  Sci  Dec;  952:23-35 


Abstracts: 
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•  AACR  Meeting,  San  Francisco,  April,  2000 

•  International  Congress  of  Endocrinology  Meeting,  Sydney,  Australia,  October,  2000 

•  Hormones  and  Cancer  Meeting,  Port  Douglas,  Australia,  November  2000 
Presentations: 

•  AACR  Meeting,  San  Francisco,  April,  2000.  Oral  presentation 

•  International  Congress  of  Endocrinology  Meeting,  Sydney,  Australia,  October,  2000.  Oral 
presentation 

•  Hormones  and  Cancer  Meeting,  Port  Douglas,  Australia,  November  2000.  Poster. 

•  Breast  Disease  Research  Group,  Baylor  College  of  Medicine,  March  20,  2002.  Seminar. 
Informatics: 

•  Full-length  G.B7  sequence  submitted  to  Genbank  database,  under  the  accession  #  AY035343 
Funding  applied  for,  based  on  work  supported  bv  this  award: 

We  have  used  the  work  presented  in  this  report  to  successfully  apply  for  a  renewal  of  NIH  Breast 
Cancer  Program  Project  grant  CA  64255,  which  will  commence  in  May  2002  under  grant  # 
CA64255-06A 

Summary  and  Conclusions 

The  E+P-treated  rat  is  a  well-established  model  to  study  the  effects  of  parity-induced  protection 
against  mammary  cancer.  In  spite  of  this  fact,  little  is  known  of  the  specific  mechanisms 
governing  pregnancy-specific  developmental  changes  in  the  mammary  gland.  We  have 
hypothesized  that  the  normal  hormonal  milieu  that  is  present  during  pregnancy  results  in 
persistent  changes  in  the  molecular  pathways  governing  cell  fate  in  a  defined  population  of  cells 
in  the  mammary  gland.  These  changes  accordingly  dictate  the  type  of  response  that  is  elicited  by 
subsequent  exposure  to  hormones  or  chemical  carcinogens.  A  critical  aspect  in  understanding 
these  processes  is  the  elucidation  of  target  genes  for  E  and  P  in  the  mammary  gland.  Such 
information  is  necessary  for  determining  how  the  molecular  pathways  involved  in  normal 
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mammary  development  and  tumorigenesis  converge  with  systemic  hormones  to  mediate  parity- 
specific  protection. 

Previous  studies,  utilizing  a  rational  approach  to  study  known  targets  associated  with 
proliferation  and  differentiation  in  the  mammary  gland,  have  met  with  limited  success  in 
revealing  the  molecular  pathways  involved  in  conferring  the  refractory  state.  Therefore,  we 
employed  an  alternative  strategy  to  identify  novel  targets  for  E  and  P  action  in  the  rat  mammary 
gland.  Using  this  Wistar-Furth  rat  model,  in  conjunction  with  subtractive  suppressive 
hybridization  (SSH)  methodologies,  we  have  identified  several  markers  that  are  persistently  up- 
regulated  in  response  to  prior  exposure  to  E  and  P  in  the  mammary  gland.  The  expression  of  two 
of  these  markers,  RbAp46  and  a  novel  gene  that  specifies  a  noncoding  RNA  (G.B7)  has  been 
characterized  further  and  the  results  of  this  work  is  described  in  a  recent  publication  (14).  We 
have  also  developed  a  real-time  PCR  assay  to  examine  the  expression  of  candidate  genes  at 
different  stages  of  development  and  in  response  to  differentiation  or  protection  of  the  gland. 

In  addition,  we  are  developing  methods  for  culturing  mammary  epithelial  cells  from  rats 
and  transducing  them  with  a  lentivirus  vector  expressing  a  gene  of  interest.  Transduced  cells 
will  then  be  transplanted  into  the  cleared  fat  pads  of  virgin  rats,  to  produce  reconstituted 
mammary  glands,  over-expressing  a  gene  of  interest.  These  experiments  are  essential  to 
advancing  our  understanding  of  the  role  of  these  markers  in  normal  mammary  gland 
development  and  in  response  to  carcinogen  challenge. 

What  are  the  implications  of  this  work? 

The  observation  that  pregnancy  confers  protection  against  breast  cancer  has  lead  to  the  proposal 
that  breast  cancer  may  be  prevented  by  pharmacologically  mimicking  the  specific  endocrine 
events  at  a  defined  window  of  development.  However,  the  advancement  of  such  technology  is 
currently  limited  by  our  ignorance  of  the  complex  molecular  and  cellular  events  that  are 
responsible  for  this  developmentally  specific  protection.  Another  factor  impeding  the 
implementation  of  such  strategies  is  the  prolonged  nature  of  any  study  that  must  rely  on  a 
defined  clinical  endpoint  (such  as  tumorigenesis)  as  a  measure  of  the  efficacy  of  any 
chemopreventative  treatment.  In  this  respect,  the  identification  of  relevant  intermediate 
biomarkers  for  a  defined  population  of  resistant  and  susceptible  cells  is  of  critical  importance, 
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and  will  greatly  facilitate  the  short-term  evaluation  of  any  chemopreventative  measures  and 
could  ultimately  improve  imaging  techniques  for  assessing  breast  cancer  risk  in  women. 

The  elucidation  of  markers  that  show  persistent  changes  in  gene  expression  in  response  to 
exposure  to  E  and  P  is  critical  for  understanding  the  molecular  pathways  that  are  altered  in  the 
parous  gland  and  modulate  the  response  of  the  gland  to  further  proliferative  stimuli.  In  this 
study,  we  have  identified  a  number  of  such  markers  that  warrant  further  study.  These  results 
provide  the  first  support  on  the  molecular  level  for  the  hypothesis  that  hormone-induced 
persistent  changes  in  gene  expression  are  present  in  the  involuted  mammary  gland  and  can  affect 
the  response  of  mammary  epithelial  cells  to  carcinogenic  insults. 

In  addition  to  the  identification  of  suitable  markers  for  following  these  parity-specific 
changes  in  the  mammary  gland,  it  is  critical  to  develop  a  method  of  overexpressing  these  genes 
in  the  rat  mammary  gland  to  enable  us  the  study  the  effects  of  their  overexpression  in  a  parity 
independent  manner.  To  this  end,  we  are  developing  a  strategy  for  transducing  mammmary 
epithelial  cells  from  rats  and  producing  reconstituted  mammary  glands  expressing  a  gene  of 
interest.  This  approach  will  enable  us  to  examine  the  role  of  a  single  molecule  (such  as 
RbAp46I.  in  conferring  protection  against  carcinogenesis,  in  the  absence  of  the  altered  hormonal 
status  that  is  present  during  pregnancy  and  without  germline  manipulation  of  the  host  animal. 
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Epidemiological  studies  have  consistently  shown 
that  an  early  full-term  pregnancy  is  protective 
against  breast  cancer.  We  hypothesize  that  the 
hormonal  milieu  that  is  present  during  pregnancy 
results  in  persistent  changes  in  the  pattern  of  gene 
expression  in  the  mammary  gland,  leading  to  per¬ 
manent  changes  in  cell  fate  that  determine  the 
subsequent  proliferative  response  of  the  gland.  To 
investigate  this  hypothesis,  we  have  used  suppres¬ 
sion  subtractive  hybridization  to  identify  genes  that 
are  persistently  up-regulated  in  the  glands  of  E- 
and  progesterone  (P)-treated  Wistar-Furth  rats 
28  d  after  steroid  hormone  treatment  compared 
with  age-matched  virgins.  Using  this  approach,  a 
number  of  genes  displaying  persistent  altered  ex¬ 
pression  in  response  to  previous  treatment  with  E 
and  P  were  identified.  Two  markers  have  been 


characterized  in  greater  detail:  RbAp46  and  a  novel 
gene  that  specifies  a  noncoding  RNA  (designated 
G.B7).  Both  were  persistently  up-regulated  in  the 
lobules  of  the  regressed  gland  and  required  previ¬ 
ous  treatment  with  both  E  and  P  for  maximal  per¬ 
sistent  expression.  RbAp46  has  been  implicated  in 
a  number  of  complexes  involving  chromatin  re¬ 
modeling,  suggesting  a  mechanism  whereby  epi¬ 
genetic  factors  responsible  for  persistent  changes 
in  gene  expression  may  be  related  to  the  determi¬ 
nation  of  cell  fate.  These  results  provide  the  first 
support  at  the  molecular  level  for  the  hypothesis 
that  hormone-induced  persistent  changes  in  gene 
expression  are  present  in  the  involuted  mammary 
gland.  ( Molecular  Endocrinology  15:  1993-2009, 
2001) 


BREAST  CANCER  IS  one  of  the  most  frequently 
diagnosed  cancers  affecting  women  in  the  West¬ 
ern  world.  Despite  advances  in  the  technologies  for 
the  diagnosis  and  treatment  of  breast  cancer  (1),  our 
understanding  of  the  etiological  factors  contributing  to 
the  development  of  this  disease  is  limited  (2).  How¬ 
ever,  one  fundamental  concept  arising  from  studies  of 
mammary  gland  biology  is  the  understanding  that  the 
endocrine  processes  involved  in  normal  mammary  de¬ 
velopment  and  carcinogenesis  are  intrinsically  related 
(3).  One  of  the  most  frequently  cited  examples  of  this 
principle  is  the  role  of  reproductive  history  and  breast 
cancer  risk.  There  is  significant  evidence  that  the  tim¬ 
ing  of  normal  developmental  events  such  as  men- 
arche,  menopause,  and  parity  and  the  duration  of  lac¬ 
tation  have  a  significant  impact  on  an  individual’s 
susceptibility  to  breast  cancer  (4, 5).  In  particular,  there 
is  strong  epidemiological  evidence  that  women  who 
experience  a  full-term  pregnancy  early  in  their  lives 
have  a  significantly  reduced  risk  of  developing  breast 
cancer  (5-7).  This  is  substantiated  by  studies  in  sev¬ 
eral  rodent  models  that  demonstrate  that  early  preg¬ 
nancy  confers  resistance  to  tumorigenesis  in  animals 
treated  with  chemical  carcinogens  (8-1 1).  These  stud- 

Abbreviations:  AMV,  Age-matched  virgin;  CoA,  coenzyme 
A;  P,  progesterone;  PPZ,  perphenazine;  Rb,  retinoblastoma; 
SSC,  standard  saline  citrate;  SSH,  subtractive  suppressive 
hybridization. 


ies  also  revealed  a  strong  age-related  component  to 
the  timing  of  developmental  events  and  exposure  to 
carcinogens  (12-14).  Thus,  the  regulation  of  normal 
developmental  events  governing  the  timing  and  period 
of  exposure  to  ovarian  hormones  appears  to  be  a 
critical  factor  in  determining  the  risk  of  developing 
breast  cancer.  Furthermore,  there  appears  to  be  a 
functional  interplay  between  the  endocrine  processes 
controlling  mammary  development  and  carcinogenesis. 

The  observation  that  the  protective  effect  of  an  early 
full-term  pregnancy  can  be  accurately  reproduced  in 
rodents  has  led  to  the  development  of  defined  animal 
models  for  studying  this  parity-related  phenomenon. 
Huggins  etal.  (1 5)  initially  demonstrated  that  treatment 
with  the  ovarian  hormones  E  and  progesterone  (P) 
could  inhibit  tumorigenesis  in  rats  after  previous  ex¬ 
posure  to  chemical  carcinogens.  Since  then,  numer¬ 
ous  investigators  have  extended  these  observations  to 
show  that  hormonal  manipulation  (by  treatment  with  E 
and  P  or  human  CG),  either  before  or  immediately  after 
carcinogen  exposure,  can  inhibit  carcinogenesis  by 
inducing  a  refractory  state.  These  studies  have  been 
performed  in  a  variety  of  rat  and  mouse  strains  (11, 
16-23).  Together,  these  findings  provide  strong  sup¬ 
port  for  the  utility  of  the  rodent  model  to  define  the 
molecular  mechanisms  by  which  a  defined  hormonal 
regimen  can  mimic  the  protective  effect  of  pregnancy. 

Despite  the  wealth  of  literature  supporting  the  role  of 
endocrine-mediated  processes  in  parity-related  re- 
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fractoriness,  little  is  known  about  the  molecular  mech¬ 
anisms  governing  pregnancy-specific  developmental 
changes  in  the  mammary  gland.  Mammary  gland  de¬ 
velopment  is  mediated  by  a  complex  interaction  be¬ 
tween  systemic  hormones  and  local  growth  factors 
(24),  which  is  in  turn  modulated  by  the  topography  of 
the  cells  receiving  the  stimuli  (2).  E  and  P  appear  to  be 
key  players  in  these  processes  (25).  However,  al¬ 
though  distinct  morphogenic  functions  have  been  as¬ 
sociated  with  these  hormones,  the  molecular  path¬ 
ways  that  are  elicited  in  response  to  the  combined 
effect  of  E  and  P  signaling  remain  to  be  elucidated  fully 
(24).  Furthermore,  it  seems  likely  that  the  signaling 
pathways  induced  by  these  hormones  vary  depending 
on  the  context  of  the  population  of  cells  receiving  the 
stimuli.  Therefore,  we  hypothesize  that  the  normal  hor¬ 
monal  milieu  that  is  present  during  pregnancy  results 
in  persistent  changes  in  the  molecular  pathways  gov¬ 
erning  cell  fate  in  a  defined  population  of  cells  in  the 
mammary  gland.  Accordingly,  these  changes  dictate 
the  type  of  response  that  is  elicited  by  subsequent 
exposure  to  hormones  or  chemical  carcinogens.  A 
critical  aspect  in  understanding  these  processes  is  the 
elucidation  of  target  genes  for  E  and  P  in  the  mammary 
gland.  Such  information  is  necessary  to  determine 
how  the  molecular  pathways  involved  in  normal  mam¬ 
mary  development  and  tumorigenesis  converge 
with  systemic  hormones  to  mediate  parity-specific 
protection. 

Previous  studies,  using  a  rational  approach  to  study 
known  targets  associated  with  proliferation  and  differ¬ 
entiation  in  the  mammary  gland,  have  met  with  limited 
success  in  revealing  the  molecular  pathways  involved 
in  conferring  the  refractory  state  (11,  21).  Therefore, 
we  used  an  alternative  strategy  to  identify  novel  tar¬ 
gets  for  E  and  P  action  in  the  rat  mammary  gland. 
Conditions  for  the  hormonal  regimen  were  based  on 
previous  studies,  which  defined  the  minimal  doses  of 
E  and  P  required  to  induce  a  level  of  morphological 
differentiation  equivalent  to  that  observed  during  a 
full-term  pregnancy  (11).  These  studies  have  shown 
that  previous  treatment  with  E  and  P  confers  resis¬ 
tance  to  chemical  carcinogen- induced  tumorigenesis 
in  Wistar-Furth  rats  in  a  reproducible  and  statistically 
significant  manner.  By  using  an  experimental  para¬ 
digm  involving  the  administration  of  E  and  P,  we  have 
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been  able  to  circumvent  the  difficulties  that  are  likely  to 
be  encountered  by  attempting  to  induce  synchronized 
pregnancies  within  an  age-matched  cohort  of  animals. 
Furthermore,  this  approach  is  likely  to  give  a  clearer 
picture  of  the  molecular  pathways  that  are  influenced 
by  E  and  P  during  normal  development  and  tumori¬ 
genesis,  The  45-d-old  rat  was  used  as  model  for  the 
initiation  of  hormonal  treatment  because  it  represents 
a  stage  of  development  that  is  analogous  to  that  of 
humans  both  in  the  rapid  development  of  the  mam¬ 
mary  epithelium  during  puberty  and  in  the  maximum 
susceptibility  to  carcinogenesis  (25).  Using  this  Wistar- 
Furth  rat  model,  in  conjunction  with  subtractive  sup¬ 
pressive  hybridization  (SSH)  methodologies,  we  have 
identified  several  markers  that  are  persistently  up- 
regulated  in  response  to  previous  exposure  to  E  and  P 
in  the  mammary  gland.  The  expression  of  two  of  these 
markers,  RbAp46  and  a  novel  gene  that  specifies  a 
noncoding  RNA  (G.B7),  has  been  characterized  in  fur¬ 
ther  detail. 


RESULTS 

Morphological  Effects  of  Hormonal  Treatments 

To  simulate  the  protective  effect  of  a  full-term  preg¬ 
nancy,  45-d-okJ  Wistar-Furth  rats  were  treated  with  a 
defined  dose  of  E  and  P  as  described  in  Materials  and 
Methods  (Fig.  1  A).  In  the  expression  studies  presented 
below,  we  used  perphenazine  (PPZ)  treatment  to  dis¬ 
tinguish  between  the  processes  involved  in  differenti¬ 
ation  and  protection  attributable  to  E+P  or  pregnancy. 
Using  the  experimental  paradigm  described  previ¬ 
ously  (22),  45-d-old  Wistar-Furth  rats  were  treated 
with  PPZ  for  a  period  of  3  wk.  This  treatment  results  in 
an  increased  serum  level  of  PRL  and  P,  but  not  E, 
which  permits  differentiation  but  does  not  confer  pro¬ 
tection  against  carcinogenesis  (22).  Hormonal  stimu¬ 
lation  (by  either  E+P  or  PPZ  treatment)  was  verified  by 
whole-mount  analysis  of  mammary  glands  at  the  con¬ 
clusion  of  the  21 -d  treatment  period.  In  both  cases, 
mammary  glands  from  6-,  12-,  and  18-d  pregnant 
animals  were  used  as  positive  controls  for  pregnancy- 
specific  differentiation  of  the  gland.  Figure  IB  shows 
representative  whole-mount  analyses  of  mammary 


Fig.  1 .  Mammary  Gland  Morphology  After  Different  Hormonal  Paradigms 
A,  Experimental  paradigm  for  studying  the  protective  effect  of  E  and  P  treatment.  Forty-two-day-old  Wistar-Furth  rats  were 
administered  E2  benzoate  (E2B)  by  sc  injection  and  then  a  sc  pellet  of  E+P  3  d  later.  The  pellets  were  replaced  after  10  d  to 
provide  hormone  treatment  for  a  total  of  21  d.  AMV  animals  were  given  blank  pellets.  After  21  d,  the  pellet  was  removed  and  the 
mammary  gland  was  allowed  to  return  to  a  resting  stage,  approximating  a  28-d  cycle  of  involution.  At  the  end  of  the  28-d 
involution  period,  mammary  glands  were  removed  for  investigation.  B,  Morphological  analysis  after  hormonal  stimulation. 
Whole-mount  preparations  of  inguinal  number  4  mammary  glands  of  66-d-old  rats  after  21  d  of  treatment  with  E+P  (i)  or  PPZ  (H). 
A  mammary  gland  from  a  1 2-d  pregnant  animal  was  used  as  a  positive  control  (iii)  to  assess  to  level  of  morphological  development 
achieved  by  the  hormonal  regimen.  A  66-d-old  AMV  animal  was  used  as  a  negative  control  for  mammary  development  in  the 
absence  of  exogenous  hormones  (iv).  Arrows  indicate  morphological  differences  between  the  E+P  and  PPZ  treatments.  C, 
Morphological  analysis  of  mammary  glands  after  hormonal  stimulation  and  involution.  Whole-mount  preparations  of  number  4 
mammary  glands  from  95-d-old  animals  treated  with  either  E+P  (i)  or  PPZ  (iii)  for  21  d  and  after  a  28-d  cycle  of  involution.  AMV 
animals  received  either  a  blank  pellet  containing  no  hormones  (ii)  or  vehicle  alone  (as  a  control  for  the  PPZ  treatment;  iv). 
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glands  from  12-d  pregnant,  21  -d  E+P-  and  PPZ- 
treated  rats,  and  age-matched  virgin  (AMV)  rats,  in  the 
E+P-  and  PPZ-treated  gland  (Fig.  IB,  i  and  il),  we 
observed  a  level  of  differentiation  that  was  comparable 
to  that  of  a  12-d  pregnant  animal  (Fig.  IB,  iii).  At  the 
gross  structural  level,  these  glands  appear  to  be  very 
similar,  although  a  few  minor  differences  were  ob¬ 
served.  In  particular,  the  E+P-treated  mammary  gland 
(Fig.  1 B,  0  appeared  to  exhibit  a  slightly  greater  degree 
of  secondary  and  tertiary  branching  than  the  PPZ- 
treated  animals  (Fig.  I  B,  b)  but  slightly  reduced  lobu- 
loaiveolar  development  compared  with  the  mammary 
glands  isolated  at  d  12  of  pregnancy  (Fig.  IB,  iii)  or 
from  the  PPZ-treated  animal  (Fig.  IB,  ii).  The  mam¬ 
mary  glands  from  21 -d  PPZ-treated  rats  exhibited  a 
level  of  lobuloalveolar  development  that  was  similar  to 
the  mammary  glands  from  12-d  pregnant  animals  (Fig. 
IB,  ii  and  iii,  respectively)  although  with  less  side 
branching  than  the  glands  from  12-d  pregnant  ani¬ 
mals.  In  addition.  21  -d  E+P-stimulated  rats  expressed 
/3-casein  at  a  level  that  was  readily  detectable  by 
Northern  analysis  using  total  RNA  (data  not  shown).  In 
contrast,  mammary  glands  from  AMV  control  animals 
receiving  blank  pellets  showed  no  significant  alveolar 
differentiation,  suggesting  that  alveolar  development 
in  the  E+P-treated  and  PPZ-treated  groups  was  a 
result  of  hormonal  stimulation.  Previous  studies  dem¬ 
onstrated  that  treatment  with  either  E+P  or  PPZ  con¬ 
fers  a  high  level  of  proliferation  and  lobuloalveolar 
differentiation  to  the  mammary  glands  of  treated  ani¬ 
mals  (22).  However,  hormonal  regulation  has  not  been 
examined  beyond  the  gross  morphological  level  to 
determine  whether  the  level  of  development  is  the 
same  in  both  cases. 

At  the  end  of  the  21  -d  treatment  period,  hormonal 
stimuli  were  withdrawn  and  the  remaining  animals 
were  allowed  to  undergo  involution  for  28  d.  This  rest¬ 
ing  phase  mimics  the  period  of  involution  that  occurs 
after  a  normal  pregnancy  and  lactation  and  represents 
a  period  of  extensive  tissue  remodeling.  At  the  end  of 
involution,  the  gland  generally  reaches  a  quiescent 
state,  resembling  that  of  the  mature  virgin  animal  (2). 

.  Figure  1C  shows  whole- mount  analysis  of  28-d  invo¬ 
luted  mammary  glands  from  95-d-old  animals  after 
previous  treatment  with  either  E+P  (i)  or  PPZ  (iii)  and 
their  age-matched  controls  (ii  and  iv,  respectively).  As 
reported  in  previous  studies  (11),  similar  morphologi¬ 
cal  characteristics  were  observed  in  all  of  the  treat¬ 
ment  groups,  indicating  that  there  were  no  persistent 
structural  differences  as  a  consequence  of  hormonal 
stimulation  of  the  gland. 

Differentially  Expressed  Genes  Identified  by  SSH 

To  help  elucidate  the  mechanisms  governing  this  hor¬ 
mone-dependent,  parity-related  phenomenon,  we  ini¬ 
tially  screened  for  molecular  markers  of  an  altered  cell 
population  or  for  changes  in  the  subset  of  genes  that 
were  persistently  expressed  as  a  result  of  previous 
hormonal  stimulation.  Using  the  experimental  para¬ 


digm  described  above,  poly(A+)  RNA  was  isolated 
from  pooled  total  RNA  samples  corresponding  to  ei¬ 
ther  E+P-treated  or  AMV  control  glands  after  a  28-d 
involution  period.  SSH  was  then  used  to  generate  an 
E+P-subtracted  library  (as  described  in  Materials  and 
Methods).  The  library  was  screened  by  differential  re¬ 
verse  Northern  blot  analysis  (26)  to  identify  genes  that 
are  differentially  induced  by  this  hormone  treatment. 
An  example  of  the  screening  method  is  presented  in 
Fig.  2,  illustrating  the  utility  of  this  approach  for  iden¬ 
tifying  differentially  expressed  markers  by  a  difference 
in  hybridization  signal  intensity  on  duplicate  filters. 
Figure  2A  shows  one  set  of  96  clones,  amplified  by 
PGR,  arrayed  by  high-density  gel  electrophoresis,  and 
photographed  to  ensure  even  loading  between  dupli¬ 
cates.  Figure  2,  B  and  C,  shows  duplicate  membranes, 
each  containing  DNA  corresponding  to  the  complete 
set  of  markers  shown  in  Fig.  2 A.  Duplicate  membranes 
were  screened  by  reverse  Northern  analysis  using 
probes  corresponding  to  either  the  E+P-subtracted 
(Fig.  2B)  or  virgin-subtracted  (Fig.  2C)  amplified  cDNA 
used  in  the  construction  of  the  SSH  library  (as  de¬ 
scribed  in  Materials  and  Methods).  This  procedure  re¬ 
sulted  in  the  identification  of  several  markers  that  were 
characterized  in  greater  detail  in  Fig.  3.  These  markers 
are  indicated  by  boxes  and  include  both  low-abun¬ 
dance  genes,  such  as  RbAp46,  Prl-1,  PP18,  and 
hnRNP  At  (Fig.  2,  B  and  G,i,  ii,  iii,  and  v,  respectively), 
as  well  as  more  highly  expressed  markers,  such  as 
transferrin  and  «~casein  (Fig.  2,  B  and  C,  iv  and  vi, 
respectively). 

Despite  the  morphological  similarity  of  the  tissues 
used  in  this  procedure,  24%  of  the  864  randomly 
selected  amplified  markers  appeared  to  be  differen¬ 
tially  expressed,  based  on  the  results  of  the  high- 
density  differential  reverse  Northern  blot  analysis  (Ta¬ 
ble  1).  Selected  differentially  expressed  Clones  (n  = 
203)  were  sequenced,  and  their  identities  were  deter¬ 
mined  using  the  BLAST  homology  search  algorithms 
and  the  public  sequence  databases  available  through 
the  National  Center  for  Biotechnology  Information. 
Clones  were  classified  as  known,  homologous,  un¬ 
known,  or  novel  depending  on  the  stringency  of  the 
search  criteria  and  the  databases  from  which  the  se¬ 
quence  matches  were  found.  A  summary  of  this  anal¬ 
ysis  is  presented  in  Table  1 .  Thus,  according  to  the 
criteria  detailed  in  the  footnotes  to  Table  1,  61.5%  of 
the  markers  sequenced  corresponded  to  known 
genes,  10.8%  corresponded  to  homologous  genes, 
and  27.5%  corresponded  to  genes  of  unknown 
function. 

A  summary  of  50  known  genes  identified  using  this 
strategy  is  presented  in  Table  2.  These  genes  fall  into 
several  distinct  categories:  those  involved  in  normal 
metabolism  and  homeostasis  of  the  gland  (metabolic 
enzymes,  transport  molecules);  those  involved  in  cell¬ 
cell  contact  and  the  extracellular  matrix;  and  regula¬ 
tory  factors  (signaling  molecules,  transcription  factors, 
etc.).  A  number  of  these  genes  (such  as  RbAp46, 
Napl ,  and  CDC42)  encode  molecules  of  known  rele- 
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Fig.  2.  Differential  Screening  of  the  E+P  SSH  Library 

An  illustration  of  the  screening  procedure  used  to  select  markers  that  are  up-regulated  in  response  to  previous  exposure  to 
E+P.  Amplified  inserts  from  96  randomly  selected  clones  were  arrayed  in  duplicate  using  high-density  gel  electrophoresis  (A). 
Gels  were  stained  with  ethidium  bromide  to  ensure  equal  loading  and  then  transferred  onto  a  charged  nylon  membrane.  Duplicate 
filters  were  screened  using  differential  reverse  Northern  analysis,  in  which  pPJdATP-labeled  probes  of  equal  specific  activity, 
generated  from  either  the  E+P  subtracted  (virgin  driver)  cDNA  (B)  or  virgin  subtracted  (E+P  driver)  cDNA  (C),  was  hybridized  under 
stringent  conditions.  The  resulting  signal  was  detected  by  exposing  the  filters  to  BioMAX  MR  film  for  6  h  to  2  d.  Six  differentially 
expressed  markers,  identified  using  this  method,  are  highlighted  by  boxes:  i)  RbAp46;  ii)  Prl-1;  iii)  PP1S;  iv)  transferrin;  v)  hnRNP 
A1;  vi)  «-casein.  The  relative  positions  of  these  cDNAs,  in  the  absence  of  a  hybridization  signal,  are  indicated  in  Fig.  2C. 


vance  to  cellular  proliferation  pathways  (27-29).  Sur¬ 
prisingly,  several  well  characterized  markers  of  mam¬ 
mary  differentiation  (such  as  the  milk  protein  genes, 
transferrin,  a-casein,  and  K-casein)  were  also  identi¬ 
fied  using  this  screen  and  displayed  persistent  expres¬ 
sion  in  E+P-treated  glands.  These  markers  were 
highly  represented  in  the  subset  of  clones  that  were 
selected  for  sequencing  and  accounted  for  18.2%, 
3.9%,  and  2.9%,  respectively,  of  the  total  number  of 
clones  sequenced.  These  observations  suggest  that 
there  is  some  degree  of  persistent  differentiation  in  the 
mammary  epithelium  after  parity  (or  treatment  with 
E+P).  However,  whether  this  differentiation  persists  as 
a  consequence  of  direct  effects  on  the  mammary 
gland  or  is  a  consequence  of  an  altered  hormonal  axis 
in  the  E+P-treated  animal  has  not  been  determined. 
Similar  findings  involving  sustained  expression  of 
these  markers  have  also  been  reported  from  studies  in 
the  parous  mouse  (3). 


Another  marker  that  was  highly  represented  was  the 
mRNA  for  the  stearyl-coenzyme  A  (CoA)  desaturase, 
an  enzyme  that  catalyzes  the  rate-limiting  step  in  the 
biosynthesis  of  monounsaturated  fatty  acids  (30). 
Stearyl-CoA  desaturase  is  expressed  in  an  isoform- 
and  temporally  specific  manner  during  development, 
depending  on  the  tissue  of  origin  (31).  The  fact  that  it 
is  so  highly  represented  among  the  subset  of  clones 
sequenced  suggests  that  there  are  also  persistent  de¬ 
velopmental  changes  in  the  mammary  stroma  in  re¬ 
sponse  to  E+P;  however,  the  implications  of  this  ob¬ 
servation  to  the  gland  as  a  whole  have  not  been 
determined. 

SSH  has  been  used  in  a  variety  of  experimental 
contexts  for  the  identification  of  differentially  ex¬ 
pressed  genes.  The  efficiency  of  the  subtraction  pro¬ 
cedure  is  validated  by  the  fact  that  several  markers 
isolated  using  this  approach  (such  as  a-casein)  have 
been  isolated  in  an  independent  study  using  quite 
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E+P  AMV  E+P  AMV  E+P  AMV 


HBP  RbAp46  Cyclophilin 


Rg.  3.  Northern  Validation  of  the  Differential  Screening  Method 
Selected  clones,  identified  by  the  SSH  library  screen  (described  in  Fig.  2),  were  used  in  Northern  blot  analysis  to  verify 
differential  expression  in  response  to  previous  exposure  to  E+P.  Northern  analysis  was  performed  with  2  /xg  of  poly(A+)  RNA  from 
either  E+P-treated  mammary  gland  or  AMV  after  28  d  of  involution.  Blots  were  hybridized  with  cDNA  inserts  prepared  from  the 
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Table  1.  Summary  of  Rat  E+P  SSH  Library  Screening 

Screening  Step 

No.  of  Clones 
Examined 

Screened  by  PCR  and  reverse 

874 

northern,  blotting 

identification  of  differentially  expressed 

213 

markers 

Sequencing 

203 

Confirm  expression  by  Northern  blot 

21 

analysis 

Percent  of  Total 

cDNA  Clones  Sequenced 

Clones 

Sequenced 

Known  genes* 

61.5 

Transferrin 

18.2 

a-Casein 

3.9 

K-Casein 

2.9 

Stearyi-CoA  desaturase 

2.S 

Homologous  genes6 

10.8 

Unknown  function0 

27.5 

ESTs 

23.6 

Novel  genes** 

•3.9 

a  Clones  classified  as  known  genes  exhibited  ^98%  identity 
with  the  corresponding  database  entry. 
b  Clones  corresponding  to  known  genes  for  which  the  iden¬ 
tity  was  not  complete  (because  the  relevant  rat  sequence  was 
not  available)  were  classified  as  homologous  if  they  demon¬ 
strated  an  expected  value  of  <1 _s. 
c  Clones  of  unknown  function  were  classified  as  such  if  they 
matched  sequences  in  the  EST  databases  but  their  function 
was  otherwise  unknown. 

^This  subset  also  included  a  number  of  novel  genes  that 
displayed  no  significant  homology  with  any  sequences  in  the 
publicly  available  databases  after  extensive  analysis  using 
the  BLAST,  TBLASTN,  and  BEAUTY  algorithms _ 


different  approaches  to  identify  genes  that  are  differ¬ 
entially  expressed  in  response  to  parity.  In  addition, 
parallel  studies  using  SSH  to  identify  markers  that  are 
down -regulated  by  E+P  have  isolated  a  quite  different 
set  of  markers  from  those  described  in  the  present 
study  (Ginger,  M.R.,  and  J.M.  Rosen,  unpublished 
data). 

Validation  by  Northern  Analysis 

To  validate  the  SSH  screen  and  reverse  Northern  anal¬ 
ysis  (described  above),  a  selected  subset  of  markers 
were  characterized  further  by  Northern  blot  analysis  to 
confirm  their  pattern  of  expression  in  the  E+P- 
preireaied,  involuted  gland  compared  with  the  AMV 
control  tissue.  A  pooled  RNA  fraction  was  used  in 


these  analyses  to  control  for  the  effects  of  variation  in 
estrous  cycle  within  anage-matched  cohort  of  animals 
(as  described  in  Materials  and  Methods).  Northern 
blots  were  hybridized  with  radiolabeled  probes  pre¬ 
pared  from  the  partial  cDNA  fragments  isolated  from 
the  SSH  library.  The  expression  of  individual  markers 
was  determined  by  quantitative  analysts  of  the  subse¬ 
quent  hybridization  signals  using  cyclophilin  as  a  nor¬ 
malizing  control.  Twenty-one  markers  were  examined 
in  this  manner,  18  of  which  were  confirmed  as  differ¬ 
entially  expressed  based  on  quantitative  Northern 
analysis  of  their  expression  in  the  E+P-treated  gland 
vs.  AMV:  Because  the  SSH  strategy  uses  an  amplifi¬ 
cation  step  to  permit  the  detection  of  low-abundance 
markers,  it  is  possible  that  some  sequences  will  be 
preferentially  amplified,  resulting  in  false-positives  re¬ 
sults.  The  observation  that  85%  of  the  markers  exam¬ 
ined  by  Northern  analysis  were  true  positives  suggests 
that  the  screening  criteria  described  above  resulted  in 
the  stringent  selection  of  differentially  expressed 
markers. 

The  results  of  the  Northern  analysis  are  summa¬ 
rized  in  Table  3  and  highlight  several  genes  that 
might  be  putative  candidates  involved  In  the  molec¬ 
ular  pathways  that  are  targets  for  E/P-mediated 
changes  associated  with  parity-related  protection  of 
the  gland.  Representative  Northern  blots  of  markers 
demonstrating  quantitative  changes  in  E-f  P-treated 
glands  are  shown  in  Fig.  3.  Several  interesting  ob¬ 
servations  accounting  for  the  differential  gene  ex¬ 
pression  were  revealed  by  these  Northern  analyses. 
Some  of  these  genes  were  expressed  as  a  single 
transcript  and  showed  a  small  quantitative  change, 
e.g.  a-casein  (Fig.  3a,  1.8-fold  difference),  Nap  1 
(Fig.  3c,  8.4-fold  change),  cdc42  (Fig.  3d.  2.1-foid 
change),  and  RbAp46  (Fig.  3q,  3.5-fold  change). 
Others  were  detected  as  multiple  transcripts,  and 
one  or  more  of  these  transcripts  were  differentially 
expressed,  e.g.  follisiati  n  -  related  protein  (Fig.  3g, 

2.1-  to  2.3-fold  change),  D.E10  (Fig.  3h,  2.9-  to 
3.4-fold  change),  hnRNP  A1  (Fig.  3k,  3.8-fold 
change).  In  the  third  case,  several  transcripts  were 
observed  and  all  of  the  transcripts  appeared  to  be 
differentially  expressed,  e.g.  G.B7  (Fig.  31,  5-  to 

7.2- fold  change).  Whether  these  transcripts  arise  as 
differently  spliced  forms  of  the.  same  gene  product 
or  closely  related  members  of  a  gene  family  remains 
to  be  determined.  Furthermore,  the  abundance  of 
these  transcripts  identified  by  SSH  and  reverse 
Northern  blotting  varied  considerably  from  highly 
abundant  mRNAs,  such  as  those  encoding  «- 


following  SSH  clones:  a,  a-casein;  b,  protein  phosphatase  1,  subunit  8;  c,  Nck-associated  protein  1;  d,  CDC42;  e,  clone  BB5, 
which  resembles  an  EST,  f,  transferrin;  g,  fullistatln  -rehied  protein;  h,  done  D.E1  D,  which  resembles  an  EST;  i,  clone  D.F2,  which 
resembles  an  EST;  j,  decorin;  k,  hnRNP  A1;  I)  G.B7;  m,  protein  phosphatase,  Prl-1 ;  n,  steary!-CoA  desaturase;  o,  ring  zinc  finger 
protein;  p,  heme-binding  protein;  q,  RbAp46.  Exposure  times  were  as  follows:  4  h  (a  and  f),  12  h  (j),  18  h  (n),  2  d  (k  and  p),  3  d 
(o),  4  d  (c,  d,  and  g),  5  d  (e  and  q),  7-8  d  (b,  h,  i,  I,  and  m).  Cyclophilin  was  used  as  a  normalizing  control  to  ensure  even  loading 
(r).  Differentially  expressed  transcripts  are  indicated  'by  snows. 
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Table  2.  Summary  of  Known  Genes  Isolated  using  SSH 


Function 


Marker 


Accession  Number 
Known  Homologous 


Metabolic 

Stearyl  CoA  desaturase 

Cysteine  dic-xygcnase 

GTP-specific  succinyl  CoA  synthetase 
2,4-Dienoyl  CoA  reductase 

J02585.1 

033431 

D00569 

AF058956 

Nucleic  acid  metabolism 

5'  Nucleotidase 

Guanine  ■cfeaminase 

J05214 

AF02S472.2 

Milk  protein 

^-Casein 

K-Casein 

0-Casein 

J0Q710.1 

K02598 

NM.017120 

Mineral  transport/metabollsm 

Transferrin 

Heme -binding  protein 

Ceruloplasmin 

D3S330 

D30035 

L33869 

Protein  synthesis 

Ribosomal  protein  LI  8 

Rj become)  protein  S20 

40-kDa  Ribosoma!  protein 

Eukaryotic  protein  translation  initiation 
factor-5 

NM.031 102.1 

X51 537.1 

D25224 

NM.020075.1 

Protein  transport 

Sec61 

U11027 

Transcription  factor 

Zinc  finger  protein,  Pzf 

Ring  zinc  finger  protein,  Rzf 

Zinc  finger  protein 

Sp3 

SHARP-2 

AF009330 

U05343.1 

AF037205.1 

U90919 

AF062567.1 

Chromatin  remodeling 

RbAp46 

AF090306 

RNA  splicing 

hnRNP  Al 

Ar  gine/serine-riob  splicing  factor  1 1 

Ml 21 56 

XMJXb  835.1 

Ceii  cycle  control 

CDC42 

Retinoblastoma 

D25233.1 

NM„009861:1 

Signal  transduction 

Krit-1 

Nap-1 

Coilibistin  \ 

D84346 

NWL023957 

AF310134.1 

Translational  repressor 

NAT-1 

U76112 

Kinase 

JAK1 

S63728.1 

Phosphatase 

Protein  phosphatase  type  1  8 

Tyrosine  phosphatase  PrI-1 

Protein  tyrosine  phosphatase-tike 
protein,  PTPLB 

D90164 

L27843 

AF1 69286 

Extracellular  matrix 

Decorin 

FoHistatiii-reiated  protein 

Heparin  sulfate  proteoglycan  core 
protein,  Sdc2 

Neura!in-1 

Z1 2298.1 

U0o364 

NM_01 3082.1 

AF305714 

Modulation  of  cell  growth 

Mast  cell'  giewth  factor 

EGF-like  protein,  S1-5/T16 

089730.1 

044725 

Cell  adhesion 

Podocalyxin 

Neural  cell  adhesion 

Celt  adhesion  molecule 

AF1 09393.1 

AF00246 

XM-0030S4.2 

Membrane  transporter 

PMP70,  ATP-binding  cassette  family 

KC1  cotransporter 

NMJ)1 2804.1 
U55815 

Celt-surface  marker 

Epithelial  glycoprotein  antigen,  EGP-314 

■fntc-rferon-  7 -receptor 

AJ001 044.1 

U68272 

Apoptosis 

E1B  19K/BcI-2 -binding  protein,  Nip3 

AF2435154 

Protease 

YMEI-like  ATP-dependent 
metal  loprotease 

NM_01 3771.1 

casein,  which  are  easily  detected  after  only  a  few 
hours-  of  exposure  (Fig,-  3a),  to  low  abundance  mR- 
NAs,  such  as  RbAp46,  which  required  more  than  1 
wk  of  exposure  (Fig.  3q).  On  the  basis  of  this  pre¬ 


liminary  screen  of  markers  up-regulated  by  treat¬ 
ment  with  E  i  Pvt wo  genes  were  selected  for  further 
characterization  in  greater  detail:  RbAp46  and  a 
novel  clone  (designated  G.B7). 
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Table  3.  Genes  Up-Regulated  by  E+P 


Marker 

Identity 

Accession  number 

Transcript* 

Size  (kb) 

Fold 

Difference* 

Known 

Homologous* 

A.F10 

Zinc  finger  protein,  Pzf 

U05343.1 

— 

B.A6 

Transferrin 

DSS380 

2,3 

2.0 

B.B5 

EST 

AI71 2974.1 

4.2 

1.9 

3.7 

— 

B.B8 

Podocafyxin 

AF1 09393.1 

1.6 

E.D2 

Sp3 

AF062567.1 

— 

B.E3 

Stearyl-CoA  desaturase 

J02585.1 

6.0 

4.4 

4.8 

— 

3.8 

— 

B.E4 

Oecorin 

Z1 2298.1 

4.2 

1.8 

3.1 

1.4 

2.0 

— 

0.9 

— 

B;F2 

n-Caseir 

J00710.1 

1.3 

1.8 

B.G10 

Zinc  finger  protein 

U90919 

— 

C.B2 

RbAp46 

AF090306 

1.9 

3.5 

C.C6 

Tyrosine  phosphatase  Prl-1 

L27843 

3.2 

— 

2,8 

6.5 

2.0 

4.7 

C.F7 

Protein  phosphatase-1 ,  6 

D90164 

2.1 

1.4 

1.5 

4.3 

C.B10 

hnRNP  A1 

Ml  21 56 

32 

3.8 

1.8 

— 

1.3 

— 

D.B4 

Ring  zinc  finger  protein 

AF037205.1 

2.2 

4.5 

D.C6 

Nap  1 

NM_023957 

5.0 

8.4 

D.F6 

Heme-binding  protein 

D30035 

1.7 

5.3 

1.1 

— 

D.F9 

EST 

BF289700.1 

2.9 

0.27 

2.2 

H  ;c 

— 

1  .o 

1.1 

2.2 

D.G2 

Follistatin-related  protein 

U06864 

3 

2.3 

2.3 

2.1 

1.6 

— 

1.3 

— 

D.E10 

EST 

BF560064.1 

4 

3.4 

2.5 

2.9 

1.7 

--- 

G.B7 

Unknown 

6.0 

7.2 

4.0 

5.6 

2.8 

7.2 

2.3 

5.0 

1.5 

6.1 

I.E1 

CDC42 

NM.009861.1 

2.2 

2.1 

a  For  ESTs,  the  accession  number  of  the  closest  homologous  sequence  is  shown. 

*  Estimated  from  -Northern  -analysis. 

c  Fold  difference  in  E+P  vs.  AMV  is  from  quantitative  analysis  of  Northern  blots  presented  in  Fig.  3  using  cyclophilin  as  a 
normalizing  control.  _ 


Cloning  of  Full-Length  G.B7 

Clone  G.B7  corresponds  to  a  764-bp  fragment  iso¬ 
lated  in  the  SSH  screen  described  above.  Analysis  of 
the  764-nucleotide  of  sequence  available  for  this 
cDNA  revealed  no  significant  homology  with  any  se¬ 
quences  in  the  publicly  accessible  nucleotide  data¬ 
bases;  hence,  it  represented  a  potentially  novel  gene 
product.  To  further  elucidate  the  function  of  this  cDNA, 


we  constructed  and  screened  an  E+P-treated  mam¬ 
mary  gland  cDNA  library  and  isolated  several  full- 
length  clones  corresponding  to  the  different  sized 
transcripts  observed  by  Northern  analysis  (Fig.  3f). 
Three  of  these  clones  (6.3,  2.4,  and  2.2  kb)  were  se¬ 
quenced  in  their  entirety.  Analysis  of  the  resultant  se¬ 
quences  suggested  that  the  different  sized  transcripts 
observed  by  Northern  blot  analysis  arose  as  a  result  of 
differential  splicing  of  a  single  gene  product.  The  full- 


2002  Mol  Endocrinol,  November  2001, 15(11): 1993-2009  Ginger  et  al.  •  Gene  Expression  in  the  Parous  Mammary  Gland 


length  sequence  of  the  largest  of  these  clones  (6.3  kb) 

was  deposited  in  GenBank  with  the.  access Idao umber 

AY035343. 

In  an  attempt  to  determine  the  function  of  this  gene, 
we  submitted  the  full-length,  6.3-kb  sequence  to  da¬ 
tabase  searching  using-  BLAST  and  TSLASTN  search 
algorithms  (with  six-frame  translation  and  BEAUTY 
postprocessing).  The  results  of  this  analysis  did  not 
reveal  significant  homology  with  any  known  gene  or 
protein  motif.  However,  we  found  that  a  short  section 
of  the  5"-  and  3r-regions  of  this  gene  (nucleotides 
1,125-1,490  and  5,920-6,203  of  the  full-length  se¬ 
quence)  exhibited  homology  with  several .  rat  and. 
mouse  expressed  sequence  tag  (ESI)  clones.  The 
best  matches  were  to  accession  numbers  BG079981 , 
BE1 19249,  and  BG079981  and  revealed  an  identity  of 
35-97%.  to  these  regions  of  approximately  300  bp  In 
the  full-length  G.B7  sequence  in  addition,  we  per¬ 
formed  homology  searching  using  the  recently  assem¬ 
bled  human  genome  sequence  database  and  found 
four  regions  of  homology  with  human  chromosome  2. 
These  homologous  regions  corresponded  to  nucleo¬ 
tides  2,334-2,419,  2,435-2,601.  2,745-2,789,  and 
5,157-5,437  of  the  full-length  G.B7  sequence,  with 
identities  of  93%;  83%,  95%,  and  80%,  respectively 
(expected  value  <  1  “°).  This  sequence  maps  to  region 
2q33  of  the  human  genome  and  spans  a  known  chro¬ 
mosomal  .breakpoint  that  is  associated-  with,  a  number 
of  human  cancers,  including  breast  adenocarcinoma 
(32).  No  open  reading  frames  or  ESTs  have  been  iden¬ 
tified  in  the  human  sequence  encompassing  the  region 
of  homology  with  G.87.  However,  this  may  be  a  con¬ 
sequence  of  incomplete  annotation  of  the  genome 
database  or  the  failure  of  the  search  paradigms  to 
detect  a  nontranslated  RNA.  Indeed  we  have  been 
unable  to  detect  any  significant  open  reading  frames 
longer  than  approximately  200  nucleotides  in  the  full- 
length  G.B7  cDNA.  In  vitro  translation  experiments, 
with  appropriate  positive  controls  run  in  parallel  to 
validate  the  assay,  also  failed  to  detect  any  translation 
product  for  the  6.3-,  2.4-,  or  2.2-kb  forms  of  this  RNA 
(data  not  shown). 

Expression  of  G.B7 

Northern  analysis  demonstrated  a  4-  to  8-fold  induc¬ 
tion  (depending  on  transcript)  of  G.87  transcripts  in 
the  glands  of  hormone- pretreated  rats  compared  with 
those  in  AMV  controls.  Multiple-tissue  Northern  anal¬ 
ysis  suggested  that  this  gene  is  also  expressed  in  the 
testes,  liver,  and  heart  of  rats  (data  not  shown);  but  its 
hormonal  regulation  in  these  tissues  has  not  been 
determined.  Because  the  mammary  gland  contains  a 
•heterogeneous  population  of  celts  (33),.  including  fibre- 
blasts,  adipocytes,  and  several  types  of  epithelial 
cells,  it  was  important  to  determine  in  which  cell  types 
this  RNA  transcript  was  expressed.  To  investigate  the 
spatial  pattern  of  G.87  expression  In  response  to 
treatment  with  E+P  and  PPZ,  we  used  in  situ  hybrid¬ 
ization  analysis  with  a  ^P-labeled  probe  generated 


from  the  764-bp  fragment  isolated  from  the  SSH  li¬ 
brary.  As  shown-  in  Fig.  4a,  G.87  mRNA'was  Taghfy 
expressed  tn  the  epithelium,  but  not  in  the  stroma,  of 
12-d  pregnant  mammary  glands.  Furthermore,  its  ex¬ 
pression  appeared  to  be  localized  to  the  differentiated 
iobuloalveolL  in  addition,  G.B7  was  persistently  ex¬ 
pressed  in  residual  alveoli  of  the  regressed  gland  after 
E+P  treatment  and  28  d  of  involution  (Fig.  4c).  By 
comparison,  only  a  very  faint  hybridization  signal  was 
detected  in  the  PPZ-treated  gland  after  28  d  of  invo¬ 
lution  (Fig.  4e).  No  signal  was  detected  in  the  95-d-old 
AMV  gland  (Fig.  4g)  or  with  the  sense  control  (data  not 
.shown).  Interesting^,  the  expression  of  G.B7  in  the 
E+P-treated  and  12-d  pregnant  animal  was  primarily 
confined  to  the  lobuloalveolar  structures  of  the  gland, 
and  not  in  the  ducts.  Furthermore,  persistent  expres¬ 
sion  of  G.87  required  previous  exposure  to  8  and  P, 
because  PPZ  treatment  alone  did  not  induce  expres¬ 
sion.  We  have  also  examined  the  expression  of  G.B7 
by  in  situ  hybridization  using  mammary  glands  from 
parous  v*s.  nenparous  animals  after  28  d  of  involution 
and  have  observed  persistent  changes  in  the  expres¬ 
sion  of  this  marker  (data  not  shown).  This  pattern  of 
expression  suggests  that  it  represents  a  potentially 
interesting  and  biologically  relevant  marker  of  parity- 
related  protection. 

Expression  of  RbAp46 

Northern  analysis  demonstrated  that  RbAp46  is  per¬ 
sistently  up-regulated  in  mammary  glands  of  E+P- 
treated  rats  compared  with  AMV  controls.  To  investi  ¬ 
gate  the  spatial  pattern  of  its  expression,  we  again 
used  in  situ  hybridization  analysis  (Fig.  5).  This  analysis 
showed  that  RhAp46  was  expressed  at  a  high  level  in 
the  tobubalveoii  at  midpregnancy  (Fig.  5,  a  and  c).  in 
addition,  RbAp46  was  persistently  expressed  in  the 
alveolar  structures  of  the  E+P-treated  gland  after  28  d 
.  of.  involution  (Fig.  5d). By  contrast,  RbAp4S  expression 
was  barely  detectable  in  the  PPZ-treated,  28-d  invo¬ 
luted  gland  (Fig.  5e)  or  the  95-d-old  AMV  control  gland 
(Fig.  5f).  RbAp46  is  also  persistently  expressed  in  the 
mammary- -glands  of  parous  animals  after- 28  d  of  in¬ 
volution  (data  not  shown).  These  results  indicate  that 
parity,  or  previous  exposure  to  both  E  and  P,  is  re¬ 
quired  for  the  maximal  persistent  expression  of 
RbAp4S. 


DISCUSSION 

The  E+P-treated  rat  is  a  well  established  model  for  the 
study  .of.  the  effects  .of  parity-induced  protection 
against  mammary  cancer.  In  spite  of  this  fact,  fittfe  is 
known  of  the  specific  mechanisms  governing  preg¬ 
nancy-specific  developmental  changes  in  the  mam¬ 
mary  gland.  Two  prevailing  models  have  been  pro¬ 
posed  to  account  for  this  phenomenon,  tn  the  first 
model  proposed  by  Russo  and  Russo  (34),  the  refrac- 
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Fig.  4.  A  Novel  Marker,  G.B7,  Is  Expressed  in  the  Lobuloalveolar  Structures  of  the  Rat  Mammary  Gland  in  Response  to  E  and  P 
In  situ  hybridization  analysis  showing  the  localization  of  G.B7  mRNA  in  the  lobuloalveoli  of  mammary  glands  from  mid 
pregnancy  (12  d  pregnant)  (a)  and  residual  alveolar  structures  of  E+P-treated  rats  (c)  after  28  d  of  involution.  As  a  comparison, 
expression  was  barely  detectable  after  treatment  with  PPZ  and  28  d  of  involution  (e)  and  in  the  95-d  AMV  glands  (g).  Hybridization 
signals  were  captured  in  the  dark  field.  4',6-Diamidino-2-phenylindole  (DAPI)  staining  of  the  same  field  is  shown  in  parallel  (b,  d, 
f,  and  h).  No  hybridization  was  observed  with  the  sense  control  (not  shown). 


tory  state  results  from  terminal  differentiation  of  the 
mammary  epithelia  during  pregnancy.  According  to 
this  hypothesis,  terminal  differentiation  rids  the  gland 
of  a  specific  population  of  susceptible  cells  (present  in 
the  terminal  end  bud  and  terminal  ducts)  that  are  con¬ 
sidered  targets  for  tumorigenesis.  Protection  is  mani¬ 
fested  as  a  change  in  cellular  kinetics,  an  alteration  of 
the  properties  associated  with  carcinogen  uptake, 
binding,  and  metabolism,  and  an  enhanced  capacity 
for  DNA  repair  (14,  25,  34-36). 

According  to  the  second  model,  the  refractory  state 
is  intrinsic  to  the  host  and  is  mediated  by  changes  in 
the  levels  of  systemic  hormones  (such  as  GH  and 
PRL),  which  may  subsequently  down-regulate  the  ex¬ 
pression  of  receptors  for  hormones  and  growth  factors 


(21).  These  models  may  not  be  mutually  exclusive,  and 
although  both  are  attractive,  neither  has  gained  suffi¬ 
cient  mechanistic  evidence  to  confirm  its  validity.  In 
addition,  comparison  of  these  hypotheses  is  ham¬ 
pered  by  the  fact  that  although  one  uses  a  pretreat¬ 
ment  model  of  hormonal  protection,  the  other  uses  a 
posttreatment  paradigm.  In  the  first  case,  differentia¬ 
tion  is  proposed  as  the  mechanism  for  the  protection, 
whereas  in  the  posttreatment  model,  an  agent  that 
causes  morphological  differentiation  of  the  gland  does 
not  confer  refractoriness  against  cancer  (22).  Further¬ 
more,  several  published  studies  have  implied  that  the 
protection  is  independent  of  the  level  of  differentiation 
achieved  by  parity  or  hormonal  stimulation  (1 0, 1 1 , 37), 
whereas  other  studies  have  contradicted  the  finding  of 
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Fig.  5.  Persistent  Expression  of  RbAp46  Requires  Previous  Exposure  to  Both  E  and  P 
In  situ  hybridization  analysis  showing  the  localization  of  RbAp46  mRN A  in  the  mammary  gland  at  1 2  d  of  pregnancy  (a  and  c) 
and  in  E+P-treated  (d)  and  PPZ-treated  (e)  rats  and  AMV  (f)  after  28  d  of  involution.  Magnification,  x4  (a,  d,  e,  and  f)  and  X10 
(c).  No  hybridization  was  observed  with  the  sense  control  (b). 


a  difference  in  cellular  kinetics  between  parous  and 
nonparous  animals  (10).  Clearly,  this  is  a  complex 
problem  that  cannot  be  fully  explained  on  the  basis  of 
existing  evidence. 

To  overcome  the  limitations  of  the  aforementioned 
models,  we  and  others  (11)  have  proposed  an  alter¬ 
native  hypothesis,  i.e.  that  the  normal  hormonal  milieu 
that  is  present  during  pregnancy  results  in  persistent 
changes  in  the  molecular  pathways  governing  cell  fate 
in  a  defined  population  of  cells  in  the  mammary  gland. 
One  possible  way  that  this  might  occur  is  through 
epigenetic  changes,  which  would  subsequently  deter¬ 
mine  the  pattern  of  gene  expression  in  a  specific  pop¬ 
ulation  of  cells  and  their  descendants.  Epigenetic 
changes,  such  as  an  alteration  in  the  methylation  sta¬ 
tus  of  promoter  sequences,  the  recruitment  of  meth¬ 
ylation  binding  proteins,  and  remodeling  of  chromatin 
structure,  would  thus  alter  the  transcriptional  profile  of 
those  cells.  Such  epigenetic  changes  frequently  ac¬ 
company  developmental  processes  and  serve  to  pro¬ 
vide  a  lasting  signal  that  restricts  the  pattern  of  gene 
expression  in  those  cells  long  after  the  inductive  event 
(such  as  hormone  or  growth  factor  signaling)  has  been 
removed  (38).  In  this  respect,  the  elucidation  of  novel 
targets  for  E  and  P  action  in  the  mammary  gland  is 
crucial  to  our  understanding  of  how  an  aspect  of  nor¬ 


mal  development  might  mediate  the  response  of  the 
organ  to  future  proliferative  signals.  We  have  selected 
the  28-d  involuted  gland  as  the  model  for  the  molec¬ 
ular  studies  presented  here  because  it  is  morpholog¬ 
ically  similar  to  that  of  mature  virgin  animals  of  the 
same  age.  Thus,  we  are  provided  with  an  appropriate 
control  to  examine  molecular  differences  in  two  tis¬ 
sues  that  are  very  similar  morphologically. 

In  the  expression  studies  presented  in  Figs.  4  and  5, 
we  used  PPZ  treatment  to  distinguish  between  the 
molecular  pathways  induced  by  either  E  and  P  or  P 
alone.  PPZ  is  a  dopamine  receptor  agonist  that  alters 
the  hormonal  axis  to  the  extent  that  both  serum  PRL 
and  P  levels  are  increased  but  E  is  not  significantly 
increased.  As  stated  above,  PPZ  caused  differentia¬ 
tion  of  the  gland,  but  not  protection,  in  a  posttreatment 
model  (22).  Whether  it  confers  protection  in  the  pre¬ 
treatment  paradigm  described  here  remains  to  be 
tested,  but  in  previous  studies  performed  in  the 
Wistar-Furth  rat  model,  both  E  and  P  were  required  for 
subsequent  protection  from  carcinogenesis  (37).  As 
we  have  shown  in  the  present  study,  both  E+P  and 
PPZ  cause  differentiation  of  the  gland.  However, 
whether  differentiation  proceeds  to  the  same  extent 
after  both  treatment  regimens  has  not  been  examined 
beyond  a  morphological  level.  The  fact  that  we  have 
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observed  slight  morphological  differences  between 
the  PP7-  and  F  4  P-treated  gland  (after  £1  d  of  stimu¬ 
lation)  suggests  that  different  molecular  pathways  are 
invoked  by  E+P-  or  PPZ-mediated  development  of 
the  gland.  Although  it  remains  controversial  whether 
differentiation  is  causal  to,  or  independent  of,  protec¬ 
tion,  treatment  with  PFZ  provides  an  important  control 
to  distinguish  between  the  signaling  pathways  in¬ 
duced  by  E  and  P  or  P  alone. 

SSH  has  been  used  in  a  variety  of  experimental 
contexts  for  the  identification  of  differentially  ex¬ 
pressed  genes.  It  has  certain  advantages  over  many 
conventional  methods,  of  gene  discovery  in  that  it  is 
capable  of  identifying  both  known  and  novel  genes  as 
well  as  low-abundance  genes.  The  last  point  is  of 
particular  interest  because  alternative  methods,  such 
as  mforoarray  analysis,  often,  fail  to  detect  Quantitative 
changes  in  low-abundance  transcripts  that  are  de¬ 
tected  by  more  sensitive  techniques  such  as  SSH  (39). 
In  many  cases,  it  is  necessary  to  use  poly(A+)  Northern 
analysis  to  verify  the  expression  of  markers  identified 
by  SSH.  In  the  present  study,  we  have  demonstrated 
the  utility  of  this  approach  for  isolating  markers  that 
are  persistently  up-regulated  by  treatment  with  E  and 
P  in  the  mammary  gland.  This  has  highlighted  a  num¬ 
ber  of  markers  of  prospective  importance  to  E-  and 
P-dependent  pathways  in  the  mammary  gland,  poten¬ 
tially  leading  to  the  determination  of  celt  fate.  Although 
some  of  these  rnarkers  represent  genes  of- unknown 
function,  others  encode  molecules  of  known  relevance 
to  the  pathways  involved  in  cellular  proliferation  and 
differentiation,  tt  is  difficult  to  reconcile  such  a  diverse 
group  of  markers  with  the  molecular  events  that  might 
be  involved  in  conferring  protection  to  the  gland.  How¬ 
ever,  it  is  possible  that  some  contribute  in  an  indepen¬ 
dent  but  stochastic  fashion  to  influence  U  se  processes 
controlling  cell  fate.  Consideration  of  their  functions 
may  lead  to  a  greater  understanding  of  the  signifi¬ 
cance  of  their  persistent  expression  in  the  E+P- 
treated  gland.  For  example,  carcinogen  treatment  of 
the  AMV  adult  gland  leads  to  a  proliferative  burst.  This 
proliferative  burst  is  attenuated  in  the  parous  and  E/P- 
treated  gland  (11).  Severn!  of  the  genes  up-regulated 
in  the  F/P-treated  gland  could  play  a  role  in  this  pro¬ 
liferative  block.  For  example,  phosphatases  play  a  crit¬ 
ical  role  in  the  regulation  of  diverse  cellular  processes, 
including  the  modulation  of  gene  expression,  celt-cycle 
progression,  and  intracellular  transport.  Folllstatin- 
related  protein  is  a  member  of  a  wider  family  of  pro¬ 
teins  that  have  been  ascribed  the  function  of  modu¬ 
lating  the  effects  of  cytokine  and  growth  factor 
signaling  (40).  In  addition,  splicing  factors  such  as 
hnRNP  Al  may  be  important  for  maintaining  the  ap¬ 
propriate  expression  of  certain  regulators  of  cell 
growth. 

We  have  characterized  two  of  these  markers  (G.B7 
and  RbAp46)  in  greater  detail  and  found  that  both  are 

persistently  expressed  in  a  specific  population  of  ep¬ 
ithelial  cells  after  treatment  with  E  !  P.  The  first  of  these 
markers  encodes  a  gene  of  as  yet  unknown  function 


that  shows  a  5-  to7-fold  induction  in  the  glands  of 
honmone-treafed  rats  .compared  with  those  in  AMV 
controls.  Toe  observation  that  it  is  homologous  with 
sequences  in  a  region  of  human  chromosome  2  and 
spans  an  area  that  is  known  as  a  chromosomal  break¬ 
point  in  a  number  of  human  cancers  is  tantalizing  and' 
merits  further  investigation.  However,  we  have  been 
unable  to  identify  a  putative  translation  product.  Be¬ 
cause  this  gene  is  clearly  expressed  as  an  mRNA 
transcript  (based  on  the  Northern  analysis  shown 
above  and  the  presence  of  related  sequences  in  the 
EST  database)  and  was  persistently  induced  after 
treatment  with  FTP,  we  speculate  that  G.87  may  in¬ 
stead  function  as  a  nontransfated  regulatory  RNA. 
There  is  some  precedence  for  this  hypothesis:  several 
studies  in  a  range  of  organisms  (including  Caenorhab - 
dills -elegans,  Drosophila,  and  mouse)  have  shown  that 
RNA  molecules  may  function  m  a  regulatory  context 
without  themselves  being  translated.  In  C.  elegans , 
two  short  noncoding  RNA  species  (fin-4  and  let-7) 
have  been  Identified  that  repress  the  function  of  sev¬ 
eral  genes  that  mediate  developmental  control  path¬ 
ways  (41).  These  RNAs  appear  to  exert  their  repressive 
effects  by  binding  to  complementary  sequences  in  the 
3 1  -untranslated  region  of  their  mRNA  targets,  thus  pre¬ 
venting  translation.  In  Drosophila,  a  family  of  noncod¬ 
ing  RNAs  encoded  by  the  roX  genes  is  involved  in 
chromatin  remodeling,  leading  to  dosage  compensa¬ 
tion-  of  the  male  X  chromosome  (42).  In  yet  another 
example,  the  RNA  molecule  SRA  acts  as  a  steroid 
receptor  coactivator  in  the  SRC-1  complex  (43).  Re¬ 
cent  studies  suggest  that  many  other  examples  of 
RNAs  with  regulatory  functions  will  be  revealed  (44). 

The  Retinoblastoma  (Rb)-associated  proteins 
(RbAp46  and  the  closely  related  molecule  RbAp48) 
constitute  part  of  a  smiail  gene  family  with  homology 
with  the  yeast  molecule  M311  (a  negative  regulator  of 
the  Ras-cAMP  signaling  pathway)  (27,  45).  Although 
originally  isolated  as  protein  components  that  bound 
to  Rb  by  affinity  chromatography  (45),  preliminary 
studies  suggest  that  these  molecules  might  play  a 
broader  function  in  the  regulation  of  such  processes 
as  cellular  proliferation  and  differentiation.  Overex¬ 
pression  of  either  RbAp46  or  RbAp48  can  substitute 
for  the  activity  of  MSI1  in  mutant  yeast  strains.  In 
addition,  RbAp46  appears  to  be  a  downstream  target 
of  the  Wilms’s  tumor  suppressor  protein  WT1V  and 
overexpression  of  this  gene  can  suppress  the  growth 
of  transfected  cells  in  culture  (46).  More  recent  studies 
have  shown  that  these  proteins  also  interact  with  the 
breast  cancer  tumor-suppressor  protein  BRCAt  (47). 
Furthermore,  RbAp46  and  RbAp43  interact  directly 
with  histones  H3  and  H4  (48)  and  are  components  of 
multisubunit  complexes  that  are  involved  in  histone 
deacetytatkm,  histone  acetylation,  nucteosome  dis¬ 
ruption,  and  nucleosome  assembly  (49-51).  In  this 
capacity,  RbAp46  appears  to  be  involved  in  both  de 
novo  acetyf  ahon/deseetytation  of  the  nascent  chroma 
tin  (possibly  leading  to  permanent  imprinting  of  spe¬ 
cific  gene  expression)  and  the  targeted  repression  of 
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gene  activation  through  its  association  with  the  Sin3/ 
HDAC  and  NuRD  complexes  (50).  Together,  these 
studies  imply  a  mode?  in  which  RbAp46,  through  its 
actions  as  an  adapter  protein,  might  serve  as  a  means 
of  recruiting  these  chromosome-remodeling  activities, 
leading  io  persistent  changes  in  gene  expression  in 
tire  glands  of  parous  animals. 

Thus,  epigenetic  changes  can  provide  an  enduring 
memory  that  predetermines  cell  fate  and  prevents  cell- 
lineage  aberrations,  leading  to  cancer  (38)-  This  offers 
a  highly  plausible  explanation  for  the  persistent 
changes  in  gene  expression  that  have  been  observed 
in  the  mammary  glands  of  parous  animals.  It  is  dear 
that  further  studies  are  necessary  to  fully  elucidate  the 
function  of  these  markers  in  the  parous  mammary 
gland.  However,  the  elucidation  of  markers  that  show 
persistent  changes  in  gene  expression  in  response  to 
exposure  to  E  and  P  Is  critical  for  under  standing  the 
molecular  pathways  that  are  altered  in  the  parous 
gland  and  modulate  the  response  of  the  gland  to  fur¬ 
ther  proliferative  stimuli.  In  this  study,  we  have  identi¬ 
fied  a  number  of  such  markers  that  warrant  further 
study.  These  results  provide  the  first  support  at  the 
molecular  level  for  the  hypothesis  that  hormone 
induced  persistent  changes  in  gene  expression  are 
present  in  the  involuted  mammary  gland  and  may 
contribute  to  the  response  of  mammary  epithelial  cells 
io  future  carcinogenic  insults. 


MATERIALS  AND  METHODS 

Animals 

Thirty- five-day-old  virgin  and  5-d  pregnant  Wistar-Furth  rats 
were  purchased  from  Harlan  Sprague  Dawfey,  fnc.  (Chicago, 
IL).  Animals  were  housed  using  approved  American  Associ¬ 
ation  for  Laboratory  Animal  Care  guidelines  in  plastic  cages 
containing  wood  chip  bedding  under  a  12-h  iight'dark  cycle 
and  were  permitted  ad  libitum  access  to  food  and  water.  AH 
experiments  were  performed  in  accordance  with  the  NIH 
guidelines  for  the  care  and  use  of  experimental  animals. 

Hormonal  Manipulation 

Forty-two-day-old  Wistar-Furth  rats  were  treated  with  a  de¬ 
fined  hormone!  mnimen  to  mimic  the  protective  effects  of  an 
early  full-term  pregnancy  using  the  experimental  paradigm 
described  previously  by  Sivaraman  et  al.  (11).  Rats  were 
divided  into  two  groups  (n  =  20  per  group):  those  receiving 
the  hormonal  stimulus  and  those  serving  as  AMV  control 
subjects.  In  both  cases,  the  rats  were  treated  with  a  priming 
dose  of  E2  benzoate  dissolved  in  sesame  oil  (2.5  pg  in  0.1  ml: 
sc)  to  synchronize  estrus  both  within  and  between  the  two 
groups.  Three  days  subsequent  to  the  E2  boost,  animals  in 
the  experimental  group  were  treated  with  E  and  P  by  sc 
implantation  of  a  peiiet  containing  20  ftg  of  £  and  20  mg  of  P 
in  a  beeswax  medium  in  the  dorsal  region  of  the  back.  Con¬ 
ditions  for  the  preparation  of  the  E/P  pellets  have  been  de¬ 
scribed  previously  (11).  AMV  control  animals  received  blank 
petrels  containing  tire  vehicle  alone.  Pefieis  were  replaced 
otter  10  d  to  provide  hormonal  stimulation  for  a  total  of  21  d. 
The  hormonal  stimulus  was  continued  for  21  d  to  mimic  the 
period  of  a  full-term  pregnancy.  After  this  phase,  the  beeswax 


pellets  were  removed  and  the  rats  were  subjected  to  a  resting 
period  of  28  d  to  allow  the  mammary  gland?  to  involute 
(Fig.  1A). 

To  study  the  effects  of  differentiation  in  the  absence  of  the 
hormonal  changes  conferred  by  pregnancy  or  E+P  treat¬ 
ment,  45-d-o!d  Wistar-Furth  rats  were  treated  with  PPZ  using 
a  modification  of  the  experimental  paradigm  presented 
above  Rats  were  divided  info  two  groups,  with  10  animals 
per  group:  those  receiving  the  hormonal  stimulus  and  those 
serving  as  AMV  control  subjects.  Rats  in  the  treatment  group 
then  received  a  sc  injection  of  PPZ  {Sigma;  5  mg/kg  in  0.03 
M  HCI)  five  times  per  wk  for  a  period  of  3  wk.  Controls 
received  the  vehicle  alone.  At  the  end  of  the  treatment  period, 
the  rats  were  again  subjected  to  a  resting  period  of  28  d  to 
allow  the  mammary  glands  to  involute. 

Tissue  Collection 

Animals  were  killed  by  an  intracardial  injection  of  0.05  ml  of 
ketamirWxyalazine/acepromazjne.  Inguinal  numb©:  4  mam¬ 
mary  glands  were  harvested  from  1 2- d  pregnant,  21-  dE-f  P 
treated,  21 -d  PPZ-treated,  28-d  involuted  (E+P,  PPZ),  and 
AMV  control  subjects.  Tissues  were  either  flash  frozen  in 
liquid  nitrogen  for  Western  and  Northern  blot  analyses  or 
fixed  in  4%  paraformaldehyde  in  PBS  at  4  C  for  1 8  h  for  In  situ 
hybridization  analysis.  For  the  12-d  pregnant  samples,  preg¬ 
nancy  was  confirmed  by  dissection  of  the  uterus  and  by  the 
presence  of  a  normal  conceptus.  For  E+P-  and  PPZ-treated 
animate,  hormonal  stimulation  of  the  gland  was  routinely 
confirmed  by  removing  the  left  number  4  mammary  gland 
immediately  after  the  21  -d  treatment  period  or  after  a  28-d 
involution  period  and  subjecting  it  to  whole-mount  analysis 
as  described  previously  {11).  In  this  case,  mammary  glands 
were  fixed  in  10%  neutral  buffered  formalin  for  24  h  and 
stained  as  described  previously  (52).  Stained  glands  were 
examined  to  ensure  that  they  displayed  morphological  de¬ 
velopment  consistent  with  the  particular  regimen. 

RNA  Isolation 

Total  RNA  was  isolated  from  frozen  tissues  by  homogeniza¬ 
tion  in  F.NAzol  B  (Tel-Test,  Friendswood,  7X)  according  to  the 
manufacturer's  instructions.  RNA  fractions  from  different  an¬ 
imals  from  the  same  group  were  combined  to  minimize  vari¬ 
ation  between  individuals  and  the  poty(A)  fraction  was  iso¬ 
lated  from  this  pooled  RNA  sample  using  the  PofyATtrsct 
mRNA  purification  system  (Promega  Ccrp.,  Madison,  Wl). 
RNA  quality  and  yield  were  determined  by  spectrophotomet- 
ric  measurement,  and  the  RNA  was  stored  at  -70  C  until 
used. 

SSH 

SSH  was  performed  using  the  POP -Set  act  oDNA  Subtraction 
Kit  (CLONTECH  Laboratories.  Inc.,  Palo  Alto,  CA)  in  accor¬ 
dance  with  the  manufacturer’s  instructions  but  with  the  fol¬ 
lowing  modifications.  Poly(A+)  RNA  was  isolated  from  pooled 
total  RNA  samples  (n  -  18  each)  corresponding  to  either 
£-f-P- treated  or  AMV  control  glands  after  a  28-d  involution 
period.  cDNA  fractions  were  synthesized  from  2  /ia  of  pofy(A) 
RNA  from  each  tissue  pool,  and  then  the  AMV  cDNA  was 
used  as  a  driver  to  subtract  molecules  common  to  both 
populations  of  RMA  from  the  E+P  tester  cDNA.  The  efficiency 
of  the  subtraction  procedure  was  monitored  by  gel  analysis 
of  the  amplified  products  before  and  after  subtraction  and  by 
depletion  of  glyceraldehyde-3-phosphate  dehydrogenase  In 
the  subtracted  vs.  the  nonsubtracted  cDNA.  Subtracted 
products  were  subjected  to  PCR-based  amplification  fet 
which  the  primary  and  secondary  PCR  conditions  were  al¬ 
tered  as  follows  to  optimize  product  formation:  in  the  primary 
PCR,  the  annealing  temperature  was  reduced  to  64  C  and  the 
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number  of  cycles  was  increased  to  30;  in  the  secondary  PCR, 
a  fin  a!  extension  cycle  of  7  mte  was  added  and  the  amplified 
products  were  “A  tailed"  by  incubation  at  72  C  for  15  min  in 
the  presence  of  0.5  U  of  Taq  polymerase  (Life  Technologies, 
lncM  Gaithersburg,  MD).  The  resultant  subtracted  amplified 
cDNA  products  were  purified  and  cloned  into  a  pGEM-T  Easy 
TA  cloning  vector  (Prcmega  Corp.).  The  ensuing  E+P  SSH 
library  was  propagated  by  transformation  into  Epicurian  Cb» 
XL2-Blue  ultracompetent  cells  (Stratagene,  La  Jolla,  CA)  ac¬ 
cording  to  the  recommendations  of  the  supplier.  Recombi¬ 
nant  E+P  SSH  cfor?es  were  selected  by  pfafihg  onto 
1 SO-mro-diameter  plates  supplemented  with  tscpropyl-l- 
thio^-D-galactopyranoside  and  5-bromo-4-ch!oro-3-indolyl- 
/3-D-galactopyranoside. 

Biffercrtriai  Screening  of  the  E+P  SSH  Library 

Colonies  (n  -  864)  from  the  E+P  SSH  library  were  selected 
at  random  and  inoculated  into  ICC  *d  of  Luria-Bcrton:  medium 
supplemented  with  amdciiHn  (100  jjug/ro J)  on  96-^weH  round- 
bottomed  microtiter  plates.  Bacteria  were  incubated  at  37  C 
overnight  on  a  shaking  platform,  and  then  3  pi  of  the  subse¬ 
quent  bacterial  culture  were  transferred  to  individual  wells  of 
a  96-wsli  PCR  plate  (Perkin-Eimer  Corp.,  Foster  City,  CA). 
Bacteria  were  lysed  by  heating  to  95  C  for  5  min,  and  the 
recombinant  DNA  inserts  were  amplified  using  nested  prim¬ 
ers  complementary  to  the  adapter  fragments  used  in  the 
It  bran/  construction,  as  described  .h  the  PCR- Select  cDNA 
Subtraction  Kit  manual  PGR  reactions  typfoaSy  contained 
0.2  mM  deoxynucleoside  triphosphate,  2.5  mw  MgCI2, 0.3  pM 
of  each  primer,  and  1  U  of  Taq  DNA  polymerase  (Life  Tech¬ 
nologies,  tec.)  in  a  total'  volume  of  20  te  the  presence  1  x 
■standard  PGR  buffer.  Taq  potemerase  was  preincubated  with 
Taqstart  antibody  (CLONTECH  Laboratories,  Inc.)  for  30  min 
at  room  temperature  before  addition  to  the  reaction  mix. 
Amplification  was  performed  using  23  cycles  each  of  94  C  for 
30  sec  and  63  C  for  3.5  min.  Equal  volumes  of  the  amplified 
products  were  arrayed,  in  duplicate,  on  high -density-format 
1 .5%  agarose  gels  (Centipede  gel  electrophoresis  chambers; 
Owl  Scientific,  Woburn,  MA),  as  described  by  von  Stein  et  al. 
(26),  and  then  transferred  onto  charged  nylon  membranes 
(Hybond  N+;  Amsrsham  Pharmacia  Biotech.  Uppsala,  Swe¬ 
den)  using  standard  protocols. 

The  resulting  duplicate  filters  were  screened  with  double- 
cONA  probes  corresponding  tc  either  reverse  or 
forward  subtracted  cDNA  prepared  during  the  library  con¬ 
struction  process.  Hybridizations  were  performed  under 
stringent  conditions  as  described  in  the  PCR-Select  Differ¬ 
ential  Screening  Kit  users  manual  (CLONTECH  Laboratories, 
.Inc,).  Hybridization  signals  were  visualized  by  exposing  the 
hybridized  filters  to  BicMAX  'MR  film  in  the  presence  of  a 
BioMAX  MS  intensifying  screen  (Eastman  Kodak  Co.,  Roch¬ 
ester,  NY),  and  the  signals  of  replicate  clones  were  com¬ 
pared.  Clones  displaying  a  differential  pattern  of  expression 
wore  selected  for  further  analysts,  and  plasmid  DMA  was 
isolated  using  a  Quantum  miniprep  kit  (Bio-Rad  Laboratories, 
Inc.,  Hercules,  CA). 

Isolation  of  FuH-Length  G.B7 

A  rat  E+P  cDNA  library  was  prepared  from  poly(A+)  RNA 

Isolated  from  26-d  Involuted  E+P-treated  mammary  glands 
and  ligated  Into  a  X7AP  Express  bacteriophage  vector  {Sirate 

agene)  in  accordance  with  the  manufacturer’s  instructions,  A 
total  of  5  x  10s  recombinant  plaques  from  the  resulting 
amplified  library  ware  screened  by  standard  protocols  using 
a  random  primed  PldCTPT  abated  cDNA  probe  corre¬ 

sponding  to  the  764-bp  G.B7  fragment  isolated  from  the  E+P 
SSH  library. 


DNA  Sequencing  and  Analysis 

Purified  plasmid  clones  were  subjected  to  dideoxy  sequenc¬ 
ing  using  an  ABI 377  automated  sequencer  (PE  Applied  Bio¬ 
systems,  Foster  City,  CA)  at  the  DNA  Sequencing  Facility  of 
the  Child  Health  Research  Center,  Baylor  College  of  Medi¬ 
cine.  Clones  from  the  SSH  library  were  sequenced  using  the 
T7  universal  primer.  Fuii-iength  G.B7  clones  were  sequenced 
using  sequence-specific  internal  oligonucleotide  primers  ob¬ 
tained  from  Integrated  DNA  Technologies  (Coralvilie,  IA).  The 
resulting  sequence  date  were  analyzed  using  the  homology 
analysis  software  (BLAST  and-  BEAUTY)  available  through 
the  National  Center  for  Biotechnology  Information  (NIH, 
Bethesda,  MD)  and  the  BCM  Search  Launcher  (Human  Ge¬ 
nome;  Sequencing  Center.  Baylor  College  of  Medidra^. 
respectively.. 

Northern  Analysis 

Tne  poiy(ATJ  RNA  fraction  from  £  P-frsafced:  .2S-d  .involuted 
and  AMV  control  mammary  glands  was  prepared  from  a 
pooled  total  RNA  sample  as  described  above.  Poly(A+)  RNA 
(2  wg/lane)  was  resolved  by  electrophoresis  on  a  1.2%  aga¬ 
rose/0.66  m  formaldehyde  gel  and  subsequently  transferred 
onto  a  charged  nylon  membrane.  The  blots  were  hybridized 
with  [a-32 P]d ATP- labeled  cDNA  probes  prepared  from  se¬ 
lected  clones  isolated  from  the  E+P  SSH  library,  stripped, 
and  reprobed  with  a  cydophliin  probe.  Hybridization  signals 
■were'  detected  by  exposing  the  fillers  let  BklMAX  -MR-  film 
(Eastman  Kodak  Co.)  in  the  presence  of  intensifying  screens. 
Several  exposure  times  were  used  to  ensure  that  the  signals 
from  indhridud  hybridizatioriS  were  in  the  linear  range  for  the 
frfrn.  films  were  scanned  using  dnnsrtemefry  (Mofeeufrr  Dy¬ 
namics,  Inc.,  Sunnyvale,  CA),  and  quantitation  was  per¬ 
formed  using  ImageQuant  1.1  software  (Molecular  Dynamics, 
Inal  The  fold  induction  of  individual  markers  was  determined 
by  normafcteg.  the  quantitative  date,  to  teat  obtained  from  the 
cyctophiJin  probe. 

in  S fttr  Hybndfzatscn 

Paraffin-embedded  sections  (7  /^m)  from  paraformalde¬ 
hyde-fixed  tissue  were  cut  and  mounted  onto  Probe-On 
Plus-charged  slides  (Fisher  Scientific.  Pittsburgh.  PA). 
Sections  were  depsraffinized,  rehydrated,  treated  with 
proteinase  K  (20  /xg/ml)  for  15  min,  postfixed  in  4%  para¬ 
formaldehyde,  and  prehybridized  for  1  h  in  hybridization 
buffer  [50%  formamide,  5x  standard  saline  citrate  (SSC; 
20  x  SSC  «  3  M  Nad.  0.3  f.i  Na3  citrate,  pH  7.0),  10% 
dexfran  sulfate.  5x  DenfvardV Vsolution;  2%  SDS4'and  TOO 
p,g/ml  denatured  salmon  sperm  DNA]  at  38  C.  [a-33P]UTP- 
Labeled  riboprobes  for  G.B7  and  RbAp46  were  transcribed 
from-  764-  and  423- bp*  fragments:,  respectively,  isolated 
f rem  the  E+P  SSH  library  cloned  mfo  pGEM-T  Easy  (Pro¬ 
mega  Corp.).  The  RbAp46  fragment  corresponded  to  nu¬ 
cleotides  1 ,279-1 ,708  of  the  fuil-length  rat  RbAp46  cDNA 
(GenBank  accession  number  AF09Q306L  Hybridization 
was  performed  at  42  C  overnight  in  the  presence  of  1  x  ID5 
cpnV/x!  radiolabeled  cRMA  probe.  Coversfips  were  re¬ 
moved  in  the  presence  of  4x  SSC  (55  C),  and  sections 
were  washed  in  2x  SSC/50%  0-mercaptoethanol  for  20 
min  at  room  temperature  and  then  digested  with  RHAss  A 
(40  p.gtevi  in  2x  SSC)  for  15  min  at  37  C,  Stringency 
washes  were  performed  in  0.1  x  SSC  for  15  min  at  42  C 
and  0.1  x  SSC  for  15  min  at  room  temperature.  Sections 
were  exposed  to  emulsion  {Eastman  Kodak  Co.)  for  3-5  -d 
and  then  mounted  in  Yeciashield  plus  4/T6-dismjdinQ-2- 
pnenylindoie  medium  (Vector  Laboratories,  Inc.,  Burlin¬ 
game,  CA). 
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